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ABSTRACT 

We present weak-lensing mass measurements of 50 X-ray luminous galaxy clusters 
at 0.15 ^ z ^ 0.3, based on uniform high quality observations with Suprime-Cam 
mounted on the 8.2-m Subaru telescope. We pay close attention to possible system¬ 
atic biases, aiming to control them at the <( 4 per cent level. The dominant source of 
systematic bias in weak-lensing measurements of the mass of individual galaxy clus¬ 
ters is contamination of background galaxy catalogues by faint cluster and foreground 
galaxies. We extend our conservative method for selecting background galaxies with 
(V — i') colours redder than the red sequence of cluster members to use a colour-cut 
that depends on cluster-centric radius. This allows us to define background galaxy 
samples that suffer ^ 1 per cent contamination, and comprise 13 galaxies per square 
areminute. Thanks to the purity of our background galaxy catalogue, the largest sys¬ 
tematic that we identify in our analysis is a shape measurement bias of 3 per cent, 
that we measure using simulations that probe weak shears upto g = 0.3. Our indi¬ 
vidual cluster mass and concentration measurements are in excellent agreement with 
predictions of the mass-concentration relation. Equally, our stacked shear profile is in 
excellent agreement with the Navarro Frenk and White profile. Our new LoCuSS mass 
measurements are consistent with the CCCP and CLASH surveys, and in tension with 
the Weighing the Giants at ~ 1 — 2a significance. Overall, the consensus at z ^ 0.3 
that is emerging from these complementary surveys represents important progress for 
cluster mass calibration, and augurs well for cluster cosmology. 
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1 INTRODUCTION 

Accurate measurements of the mass and internal struc¬ 
ture of dark matter halos that host galaxy clusters and 
groups are central to a broad range of fundamental research 
spanning cosmological parameters, the nature of dark mat¬ 
ter, the spectrum of primordial density fluctuations, test¬ 
ing gravity theory, and the formation/evolution of galax¬ 
ies and the intergalactic medium. The requirement for ac¬ 
curacy is most stringent for studies that aim to probe 
dark energy, e.g. via evolut i on of the cluster mass functio n 
(e.g. IVikhlinin et al.l 120091 : lAllen, Evrard fc Mantd 1201 ll) . 
Upcoming surveys will discover ~ 10 5 clusters and intend to 
infer the mass of the majority of these systems from scaling 
relations between mass and the observable properties of clus- 
ters le.g. IPillepich, Porciani &: Reiprichll2012l ; ISartoris et al.l 
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f E-mail: gps@star.sr.bham.ac.uk 


120151) . Notwithstanding the forecast accuracy and precision 
of other cosmological probes, the sheer number of clusters 
upon which future cosmological results will rely implies that 
per cent level control of systematic biases in the ensemble 
mass calibration of clusters will ultimately be required. 


The challenge of calibrating systematic biases in the 
ensemble galaxy cluster mass calibration at this level of 
accuracy is amplified by the fact that the normalization 
of the calibration is necessary but not sufficient for ac¬ 
curate cluster cosmology. Poorly constrained knowledge of 
the intrinsic scatter between observable mass proxies (in¬ 
cluding all “masses” measured from data) and the under¬ 
lying mass of dark matter halos that host galaxy clus¬ 
ters is a sour ce of bias in clust er-based cosmological con¬ 
straints le.g. ISmith et all 1200 3). Intrinsic scatter between 
the relevant observable properties of clusters and between 
mass measurements and underlying halo mass is there¬ 


fore a key parameter that many stud i es attempt to con- 
strain (e.g. lOkabe et al. ] l2010bl . I2014bl : iBecker fc Kravtsovl 
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201 ll; Bahe, McCarthy & Kind 20121; 

Marrone et al. 

2012 

Mahdavi et al.l 20131: Sifon et al. 2011 

: Mulrov et al. 

2014 


Rozo et al.l 20151 : Saliwanchik et alT 2015l l. Cluster mass 


measurement methods used for calibration studies must 
therefore permit measurements of individual cluster masses 
in order to characterise the full distribution of cluster mass 
around the mean relation between mass and observable 
mass proxy. Moreover, whilst stacked mass measurements 
are powerful probes of the population mean, they offer no 
useful constraints on the scatter around the mean. 

An increasing number of galaxy cluster mass calibration 


strain cluster masses 

e.g. Smith et al. 20051: Bardeau et al. 

2007 

; Okabe et al. 

2010al: Okabe. Okura & Futamase 

201C 

; (Okabe et al. 

201ll. 20131. 2014al. 20151. 20ld: 

Hoekstra et al.l 20121: 

Applegate et al. 20141; Umetsu et al. 

20141; Hoekstra et al. 201511. This is because interpreta- 


tion of the gravitational lensing signal does not require 
assumptions about the physical nature or state of the 
gravitating mass of the cluster. Therefore, despite the 
fact that individual cluster mass measurements can suffer 
appreciable biases that correlate with the observer’s viewing 
angle through asymmetric cluster mass dist ributions (e.g. 
ICorless fe Kindl2007l : iMeneehetti et al.lfeoidl . gravitational 
lensing can yield an accurate mean mass calibration of 
galaxy clusters, supported by knowledge of the scatter be- 
tween true halo mass and weak-lensing mass measuremen ts 
dBecker fe Kravtsovll201 3 ; iBahe, McCarthy fe Kingll2012l l. 

The largest samples of clusters for which weak-lensing 
observations are available are currently drawn from large- 
scale X-ray surveys and number of order 50 clusters. These 
survey s are the Local Cluster Substructure Survey (Lo- 
CuSS: lokabe et al.l 1201 .'ll : iMartino et al.l 120141 : Smith et al. 
in prep.), the Canad i an Cluster Cosmology Project (CCCP; 
iMahdavi et al.lf2013l : iHoekstra et al.ll2015h. and the Weigh¬ 
ing the Giants programme (WtG: [von der Linden et al.l 
120141 : lApplegate et al.l l2014l l. In the parlance of the Dark 
Energy Task Force, these are Stage II studies that examine 
the systematic uncertainties inh erent in using galaxy clus¬ 
ters as probes of Dark Energy dAlbrecht et al.l [20061 1. The 
LoCuSS sample is an Lx-limited sub-set of clusters from 
ROSAT All-sky Survey (RASS) at 0.15 < 2 < 0.3; the 
CCCP sample is a mixture of X-ray luminous clusters for 
which optical data are available from the CFHT archive and 
a temperature-selected sub-set of clusters from the ASCA 
survey spanning 0.15 < 2 < 0.55; the WtG sample is a 
representative flux-limited sub-set of the RASS clusters at 
0.15 < 2 < 0.7. Smaller, generally heterogeneous, samples of 
X-ray clusters are also studied, for example, by the the Clus- 
ter Lensing And Supernova Survey with Hubble (CLASH; 
IPostman 2012 • Umetsu et al.l l2014h and the 400SD 
surveys ( Isra el e t al.ll2012i l . Whilst samples of Sunyaev Zel- 
dovich (SZ) Effect detected clusters are growing rapidly, the 
weak-lensing studie s of SZ samples currently number ha nd- 
fuls of clusters (e.g. iHigh et al.ll2012l ; iGruen et al.ll20l4l . 

Currently, the target accuracy on controlling systematic 
biases in the ensemble cluster mass calibration is therefore 
set by the size of the LoCuSS, CCCP and WtG samples, 
the typical statistical measurement error of a weak-lensing 
mass measurement of an individual cluster, and the intrinsic 
scatter of weak-lensing masses around the true underlying 
halo masses. Given that our sample is not mass-selected, 


the intrinsic scatter on Mwl — Mtrue for our sample is not 
known a priori. In setting a nominal goal for control of sys¬ 
tematic biases, we therefore ignore the intrinsic scatter and 
simply adopt a typical statistical measurement error of 30 
per cent as the uncertainty on the ensemble cluster mass 
calibration that would be achieved from studying one clus¬ 
ter. This motivates a goal of ~ 30/\/50 = 4 per cent for 
control of systematic biases. Our goal in this article is to 
achieve this level of accuracy for the LoCuSS galaxy cluster 
mass calibration. Note that, by ignoring the intrinsic scatter 
in Mwl — Mtrue, this goal is more challenging than is justi¬ 
fied by the statistics of cluster mass measurement discussed 
above. 

The principal systematic biases that can affect an weak- 
lensing cluster mass measurement relate to (1) the mea¬ 
surement of faint galaxy shapes, (2) accurate placement of 
faint galaxies along the line of sight such that the sample 
of background galaxies suffers negligible contamination by 
faint cluster members and that the inferred redshift distri¬ 
bution of the background galaxies is accurate, and (3) mod¬ 
elling of the shear signal in order to infer the cluster mass. 
In the brief review of these biases that follows, a key theme 
is that the approach taken to addressing one source of bias 
can have consequences for how well other biases are con¬ 
trolled. We also briefly outline our approach to these biases 
- a unifying theme of which is to minimize the number of 
strong assumptions in our analysis. The summary that fol¬ 
lows intends to help non-experts to understand some of the 
more technical aspects of this article. 

It is common to calibrate faint galaxy shape mea- 
surement codes on t h e STEP and STE P2 simulations 
dHevmans et al.l I2OO6I : iMassev et al.l l2007l l. however the 
gravitational shear signal of clusters typically exceeds the 
shear signals injected into these simulations, and therefore 
these tests are only relevant to the cluster outskirts. For ex¬ 
ample, WtG calibrate their shape measurement code on the 
STEP2 simulations, which in part motivates them to restrict 
the range of the WtG shear profiles to projected cluster¬ 
centric radii of 0.75 — 1.5 hj 0 Mpc - i.e. attempting to avoid 
regions of the clusters at which the measur ed shear exceeds 
that i njected into the STEP simulations (lApplegate et al] 
|2014[) . IHoekstra et al.l (l2015l l recently emphasised the im¬ 
portance of carefully matching the properties of the sim¬ 
ulated data used for such tests to the observational data. 
In this article we further dev elop the shap e mea surement 
methods that we developed in lOkabe et al.l (|2013l l and ex¬ 
tend our tests of these m ethods to include realistic galaxies. 
As in lOkabe et alj J2013ll . we test our code on shear values 
upto g = 0.3, i.e. appropriate to the full range of cluster 
centric radii relevant to weak-lensing - clustercentric radii 
as small as ~ 200 h _1 kpc. 

Contamination of background galaxy samples by un- 
lensed faint cluster galaxies dilutes the measured lensing sig- 
nal and causes a systematic underestimate o f the shear (e.g. 
iBroadhurst et aPl2005l : iLimousin et al.l[2007l l. It is therefore 
of prime importance to make a secure selection of back¬ 
ground galaxies. The number density of cluster members 
is a declining function of clustercentric radius, and thus a 
radial trend in the number density of galaxies selected as 
being in the background is interpreted as evidence for con¬ 
tamination. Whilst this is qualitatively true, the quantita¬ 
tive details depend on how the gravitational magnification 
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of the cluster modifies the observed distribution of back¬ 
ground galaxies. A fter falling into disuse for a decade since 
iKneib et al.l d2003h first proposed the method, boosting the 
measured shear signal to correct statisticall y for contamina¬ 
tion has enjoyed a renaissa nce of late (e.g. lAppleeate et al.l 
120141 : lifoekstra et al.ll2015l i . This method is applied to both 
red and blue galaxies, either by excluding the red sequence 
galaxi es or simply selecting faint galaxies, as per lKneib et al.l 
J2003I 1. Due to imperfect background selection, the number 
density profile of these colour-selected galaxies is found to 
increase at small cluster-centric radii. Assuming that the 
number density profile of a pure background galaxy sam¬ 
ple is independent of radius, that is ignoring gravitational 
magnification, the lensing signal is corrected as a radial func¬ 
tion of the galaxy-count excess. This correction method is 
referred to as “boost correction”. 


We also note that photometric redshifts based on upto 
five photometric bands are becoming mor e common as a 
meth od for selecting ba ck ground ga laxies dLimous in et al.1 


20071: Gavazzi et al. 

20091: Gruen et al.l 2013i: 

ADDlegate et al. 2014 

: McClearv. dell’Antonio & Huwc 

20151: Melchior et al. |2015f). However photometric red- 


shifts based on a small number of filters are problematic 
for galaxies with blue observed colours because their 
spectral energy distribution is relatively featureless. This 
leads to the well known degeneracy between photometric 
redshifts of z < 0. 5 and z > 1.5 for blue galaxies (e.g. 
iBolzonella, Miralles fe Pellol l200ol l. This is a critical issue 
for cluster weak-le nsing studies that use blue galaxies 
(iZiparo et alJ l2015h . Furthermore, the requirement for 
photometric redshift accuracy is more stringent for cluster 
lensing studies than for most other fields, because the 
number density of cluster galaxies - that contaminate 
background galaxy samples - is a function of clustercentric 
radius. 

We have previously developed a method to select red 
background that does not assume the radial distribution of 
background galaxies and thus does not require a boost cor - 
rection to the measured shear signal (lOkabe et al.ll2013l) . 
Our method also yields a direct measurement of the frac¬ 
tion of galaxies i n the background g alaxy sample that are 
contaminants. In IZiparo et al.l (12015T ) we considered how to 
extend this method to include blue galaxies and concluded 
that additional uncertainties of including blue galaxies do 
not justify the small number of additional galaxies that we 
would gain. We therefore extend our red galaxy selection 
methods in this article, and achieve a 2.6-f old increase in 
number density of background galaxies over lOkabe et al.l - 
i.e. sufficient to measure individual cluster masses, whilst 
retaining our conservative requirement that contamination 
is not greater than 1 per cent. 

Despite the intrinsic asphericity of galaxy clusters, 
it has been shown that modelin g cluster mass distribu¬ 
tions a s spherical and following a lNavarro. Frenk fe White! 
(1 1997a| j profile yields mass meas urements that are accu¬ 
rate in the mean across a samp le dBecker fe Kravtsovll201 J : 
iBahe, McCarthy fo King|l2012l f. These results are based on 
numerical dark matter only simulations, make (well moti¬ 
vated) assumptions about the observational data available 
to an individual study, and stress the importance of fitting 
the model to the data across a well-defined radial range. 
Some observational studies implement directly the method 


desc ribed b y [Be cker &: Kravtsov! (1201 il l in their analysis 
fe.g. lApplegate et al.ll2014 1. We prefer to test our mass mod¬ 
elling scheme directly on simulations. Moreover, parameters 
that describe the shape of the density profile (generally, a 
“halo concentration parameter”) are at the same time a nui¬ 
sance parameter for the mass measurement, and a physically 
interesting parameter to extract from the data. We there¬ 
fore let concentration be a free parameter with a flat prior, 
and marginalise over concentration when measuring cluster 
mass. We also use the constraints that we derive on concen¬ 
tration to examine the mass-concentration relation. Other 
studies adopt more restrictive assumptions about halo con¬ 
centration, in part as a consequence of seeking to minimise 
contamination and shear calibration issues (see preceding 
discussion) by excluding the central cluster region from their 
analysis and modeling. 

We describe the observations and data analysis in Sec¬ 
tion [2j including photometry, shape measurements, and the 
selection of background galaxies. The mass measurements 
for individual clusters, the mass concentration relation and 
stacked lensing analysis are presented in Section [3] We dis¬ 
cuss several systematics and compare with previous weak 
lensing studies in Section [4j and summarize our conclu¬ 
sions in Section [ 0 ] We assume Ho = lOOh kms _1 Mpc _1 , 
D m o = 0.3 and Da = 0.7 through the paper. We occasion¬ 
ally use the alternative definition of the Hubble parameter 
h 70 = H o /70. 


2 DATA AND ANALYSIS 
2.1 Sample 

The sample comprises 50 clusters (Table [T]) d rawn from the 
ROSAT All Sky Survey clu ster catalogues dEbeling et al.l 
1 19981 . 120001 : iBohringer et al.l l2004lf that satisfy the crite¬ 
ria: -25° < S < +65°, n H ^ 7 x 10 2O cnW 2 , 0.15 < 2 < 
0.3, Lx/E(z) > 4.1 x 10 44 ergs^ 1 where Lx is in the 
0.1 — 2.4 keV band and E{z) = ^/(1m (1 + z) 3 + Da. The 
sample is X-ray luminosity-limited, and therefore approxi¬ 
mately mass-limited. Full details of the selection function 
are available in Smith et al. (2016, in prep.). 


2.2 Observations 

The clusters were observed with Suprime-Cam 
dMivazaki et al.l 120021 1 on the 8.2-m Subaru TelescomFTl 
on Mauna Kea. We observed in both V- and i'-bands for 
28 and 36 minutes respectively, splitting the integration 
times up into individual four minute exposures. The 
best overhead conditions were reserved for the i'-band 
observations because we use these data to measure the 
shapes of faint galaxies. The full width half maximum 
(FWHM) of point sources was routinely sub-arcsecond, 
with individual exposures often enjoying FWHM < 0.6 ,/ . 


1 Based in part on observations obtained at the Subaru Obser¬ 
vatory under the Time Exchange program operated between the 
Gemini Observatory and the Subaru Observatory. 

2 Based in part on data collected at Subaru Telescope and ob¬ 
tained from the SMOKA, which is operated by the Astronomy 
Data Center, National Astronomical Observatory of Japan. 
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Table 1. Cluster sample. 1 ) :cluster name. 2 ); cluster redshift. 3 ) :the band name for shape measurement. 4 ) :the band combinations for the 
color. 5 );the seeing size for the WL bands. 6 ) :the magnitude range for the red-band. 7 );the number density for background galaxies after 
the color selection. 8 ):the signal-to-noise ratio in the tangential shear profile. 


Name 1 ) 

Redshift 2 ) 

WL 3 ) 

color 4 ) 

Seeing 5 ) 

[arcsec] 

mag 6 ) 

[ABmag] 

^bkg”^ 
[arcmin -2 ] 

S/N 8 ) 

ABELL2697 

0.2320 

V 

i'V 

0.77 

21 - 26.0 

15.7 

6.5 

ABELL0068 

0.2546 

i' 

i'V 

0.73 

21 - 25.9 

19.0 

8.1 

ABELL2813 

0.2924 

i' 

i'V 

0.77 

21 - 26.1 

10.2 

6.7 

ABELL0115 

0.1971 

i' 

i'V 

0.71 

21 - 24.3 

7.2 

3.9 

ABELL0141 

0.2300 

i' 

i'V 

0.71 

21 - 26.2 

17.7 

5.7 

ZwC10104.4+0048 

0.2540 

V 

i'V 

0.65 

21 - 25.8 

15.5 

2.8 

ABELL0209 

0.2060 

%' 

i'V 

0.63 

21 - 24.9 

13.7 

9.4 

ABELL0267 

0.2300 

i' 

i'V 

0.61 

21 - 25.6 

20.1 

7.2 

ABELL0291 

0.1960 

i' 

i’V 

0.71 

21 - 25.9 

15.8 

5.9 

ABELL0383 

0.1883 


i'V 

0.67 

21 - 25.9 

19.8 

7.6 

ABELL0521 

0.2475 

i' 

i'V 

0.61 

21 - 24.8 

15.1 

6.2 

ABELL0586 

0.1710 

i' 

i'V 

0.83 

21 - 25.6 

6.7 

6.3 

ABELL0611 

0.2880 

i' 

i'V 

0.79 

21 - 25.8 

8.0 

5.7 

ABELL0697 

0.2820 

i' 

i'V 

0.73 

21 - 26.1 

12.8 

6.5 

ZwC10857.9+2107 

0.2347 

i' 

i'V 

0.85 

21 - 26.0 

9.2 

2.6 

ABELL0750 

0.1630 

i' 

i'V 

0.71 

21 - 25.8 

13.4 

6.2 

ABELL0773 

0.2170 

i' 

i’V 

0.57 

21 - 26.5 

21.4 

10.5 

ABELL0781 

0.2984 

V 

i'V 

0.87 

21 - 26.1 

8.4 

4.2 

ZwC10949.6+5207 

0.2140 

i' 

i'V 

0.83 

21 - 26.2 

9.9 

5.1 

ABELL0901 

0.1634 

%' 

i'V 

0.73 

21 - 26.2 

12.9 

4.0 

ABELL0907 

0.1669 

i' 

i'V 

0.73 

21 - 26.4 

13.1 

8.8 

ABELL0963 

0.2050 

/c 

IcV 

0.75 

21 - 25.7 

15.0 

7.2 

ZwCU021.0+0426 

0.2906 

i' 

i'V 

0.61 

21 - 26.3 

20.4 

6.7 

ABELL1423 

0.2130 

i' 

i'V 

0.71 

21 - 26.0 

16.5 

6.0 

ABELL1451 

0.1992 

i' 

i'V 

0.65 

21 - 26.4 

21.3 

10.2 

RXCJ1212.3-1816 

0.2690 

i' 

i'V 

0.91 

21 - 26.1 

10.2 

2.1 

ZwC11231.4+1007 

0.2290 

i' 

i'V 

0.75 

21 - 26.0 

8.7 

4.4 

ABELL1682 

0.2260 

i' 

i'V 

0.69 

21 - 26.2 

21.9 

10.8 

ABELL1689 

0.1832 

i' 

i'V 

0.87 

21 - 25.7 

8.9 

10.6 

ABELL1758N 

0.2800 

Rc 

i'B 

0.69 

21 - 25.9 

13.7 

4.8 

ABELL1763 

0.2279 

%' 

i'V 

0.77 

21 - 25.5 

11.7 

9.6 

ABELL1835 

0.2528 

i' 

i’V 

0.89 

21 - 25.2 

8.4 

7.9 

ABELL1914 

0.1712 

Rc 

i'g' 

0.61 

21 - 26.1 

11.7 

6.7 

ZwC11454.8+2233 

0.2578 

i' 

i'V 

0.81 

21 - 25.2 

6.9 

3.1 

ABELL2009 

0.1530 

i' 

i'V 

0.77 

21 - 26.1 

12.8 

6.1 

ZwCU459.4+4240 

0.2897 

Rc 

i'V 

0.57 

21 - 26.0 

12.3 

3.7 

RXCJ1504.1-0248 

0.2153 

i' 

i'V 

0.79 

21 - 24.6 

6.7 

5.2 

ABELL2111 

0.2290 

i' 

i'V 

0.89 

21 - 25.3 

8.3 

3.4 

ABELL2204 

0.1524 

i' 

i'V 

0.81 

21 - 24.9 

8.6 

7.9 

ABELL2219 

0.2281 

i' 

i'V 

0.71 

21 - 26.0 

17.2 

7.7 

RXJ1720.1+2638 

0.1640 

i' 

i'V 

0.71 

21 - 24.3 

7.5 

4.7 

ABELL2261 

0.2240 

i' 

i'V 

0.61 

21 - 26.0 

21.3 

10.8 

RXCJ2102.1-2431 

0.1880 

V 

i'V 

0.71 

21 - 25.9 

14.2 

4.0 

RXJ2129.6+0005 

0.2350 

i' 

i'V 

0.85 

21 - 25.8 

11.6 

4.6 

ABELL2390 

0.2329 

Rc 

i'V 

0.65 

21 - 26.2 

12.3 

8.9 

ABELL2485 

0.2472 

i' 

i’V 

0.67 

21 - 25.6 

17.0 

5.9 

ABELL2537 

0.2966 

V 

i'V 

0.99 

21 - 25.9 

4.9 

5.2 

ABELL2552 

0.2998 

Rc 

i'V 

0.77 

21 - 26.0 

9.9 

4.3 

ABELL2631 

0.2779 

Rc 

i'V 

0.65 

21 - 25.9 

12.4 

4.5 

ABELL2645 

0.2510 

i' 

i'V 

0.67 

21 - 26.0 

11.9 

4.8 


The 50 final stacked and reduced i'-band frames have shapes as the i'-band, and the bluer filter used for colour 

median seeing of FWHM = 0.71", with 38 of the 50 frames measurements as the F-band. 

having FWHM < 0.8" (Table [TJ. Note that we use archival 

g- and B-band dat a instead of F-ban d data for two clusters 

in common with lOkabe fc Umetsul (120081 '). Hereafter we 

refer to the redder filter in which we measure faint galaxy 
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2.3 Data Reduction 


We reduced all data using a processing pipeline based on 
the the standard reduction tasks for Sup rime-Cam, SD- 
FRRD_([Yagi_fjt ah||2002|: lOuchi et al.l 120041 . and described 
bv lOkabe et al.l 1 201031 1. The pipeline includes bias and dark 
frame subtraction, flat-fielding, instrumental distortion cor¬ 
rection, differential refraction, point spread function (PSF) 
matching, sky subtraction and stacking. The astrometric so¬ 
lution for t he final stacked frame s was calibrated relative 
to 2MASS (ISkrutskie et al.l l2006l l to sub-pixel root mean 
square (rms) precision. Photometric zero-points were cal¬ 
ibrated to stellar photometry from the Sloan Digital Sky 
Survey dEisenstein et al.l 1201 ll . SPSS), taking into accoun t 
foreground galactic extinction dSchlaflv fe Finkbeineil201l|j . 
to a rms precision of ,$ O.lmag. To cross-check the validity 
of the photometric calibration, we measured the redshift de¬ 
pendence of the colour of early-type member galaxies, that 
lie on the so-called cluster red-sequence, within 10 arcmin of 
each brightest cluster galaxy (BCG). The colour of the red 
sequence increases from (V — i ') ~ 0.8 to (V — i') ~ 1.2 
as cluster reds hift increases from z — 0.15 to z = 0.3, in 
agreement with lEisenstein et al.l (1201 il l . 


2.4 Shape measurement pipeline 

We anal yse the i'-band frames wi t h me thods intro¬ 
duced bv lKaiser, Squires fc Broadhurstj d 19951 . the “KSB-(-” 
method), using the imcat package with our modifications 
dOkabe et al.l 120131 . l2014al l. We first measure the image el- 
lipticity, e a , from the weighted quadrupole moments of the 
surface brightness of objects, and then correct the PSF 
anisotropy by solving 

eUe)=e a {e)-P^(B)q}{G), (1) 

where P a p is the smear polarizablity tensor and <j£(0) = 
(Psm)kp e *’ quantities with an asterisk denote those for stel¬ 
lar object s. The details of anis otropic PSF correction is de¬ 
scribed in lOkabe et al.l d2014al . see the Appendix). In brief, 
we selected bright, unsaturated stars in the half-light radius, 
rh, and magnitude plane to estimate the stellar anisotropy 
kernel, Note that the stars and galaxies can be clearly 
discriminated using the half-light radius. We modeled the 
variation of this kernel across sub-regions of the field of view 
by fitting second-order bi-polynom ial functions to the vec¬ 
tor 0 with iterative g-clip ping (e.g. lOkabe fe Umetsull2008l : 
lOkabe et alJl20fo3 . l2Q14ah . Although distortions at the cor¬ 
ners of the field-of-view are larger than those at the centers, 
modelling across sub-regions is sufficiently flexible to correct 
the anisotropic PSF pattern in our data. 

We tested the validity of our anisotropic PSF correction 
by measuring the auto-correlation function between stellar 
ellipticities and the cross-correlation function between stel¬ 
lar and galaxy ellipticities, before and after the correction. 
We found that (e^’ raw e^’ raw ) and (e a e^’ raw ) before the cor¬ 
rection show large positive correlations (C?(10 -4 )), and that 
(ea Tes ea res ) and (e“ r ea’ res ) after the correction are consis¬ 
tent with null correlation in individual cluster fields. Note 
that e“ r is the l.h.s. of Equation [Q In order to confirm 
how well the anisotropic PSF correction works at both the 
corners and centers, we divide the regions into the inner 
(r < 14') and outer regions (r > 14') of the fields of view 


Inner Region: r<14 [arcmin] Outer Region: r>14 [arcmin] 



• * • 

$ a = l 
$ a—2 

— fl - ** - ♦* 

— ii — a — *i- 

-jj — 

b*H 

l-f-d 


1 5 5 

9 [arcmin] 9 [arcmin] 


Figure 1. Top - The cross-correlation between raw elliptici¬ 
ties for stars an galaxies before the correction. The correlation 
at the outer region (right) is larger than that at the inner re¬ 
gion (left). Middle - The cross-correlation between residual stel¬ 
lar ellipticities and corrected galaxy ellipticities. Bottom — The 
cross-correlation between residual stellar ellipticities and reduced 
ellipticities. 


with respect to the respective BCGs (Figure [TJ. Note that 
the BCGs are located close to the center of the field of 
view in all cases. The top-left and top-right panels show the 
cross-correlation function between stellar and galaxy ellip¬ 
ticities before the correction at the inner and outer regions, 
respectively. The raw distortions at the outer region is in¬ 
deed larger than those at the inner region. The middle panel 
shows the resulting (e“ r ea res ) after the correction, which are 
consistent with null correlation both at r < 14' and r > 14'. 
The cross correlation between the residual stellar ellipticities 
and the reduced ellipticities for galaxies which are described 
in next paragraph is shown in the bottom panel. We again 
found null correlation at r < 14' and r > 14'. 

Next, we correct the isotropic smearing effect of galaxy 
shapes due to seeing and the Gaussian window function used 
for the shape measurements. The reduced ellipticity for each 
galaxy, g a , is defined by 

g a = {Pg^otpe'fj, (2) 


where (P g ) a p is the pre-seeing shear polarizability 
tensor. The measurement of ( P g ) a p is very noisy for 
individual faint galaxies bec ause of its non-linearity 
jBartelmann fe Schneideu 1200111 . The relationship between 
noise and biases in the measurement of faint galaxy 
shapes, using a variety of shape measurement algorithms, 
has been considered for both cosmic s h ear and cluster 
lensin g studies fe.g. iH irata e t al.l 12004; iKacprzak et alJ 
2012|__ Melchior fe Violal 201^__ Refregier et al.l 20121 : 


Applegate et al. 2014 : Hoekstra et alJ 2015h . The gen¬ 


eral feature of this relationship is that the bias in faint 
galaxy shape measurements typically increases as the 
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6 Okabe and Smith 


size of galaxies decreases, i.e. as signal-to-noise ratio of 
galaxies decreases. In common with s everal authors (e.g. 


Umetsu et al. 20101: 

Oguri et al.l -20121: lOkura & Futamasel 

20121: lUmetsu et al.l 

2015a|) we have found that this de- 


pendence can be reduced significantly for KSB shape 
measurement methods if the galaxies upon which the 
isotropic PSF correction is based ar e limited to thos e 
detected at high signal-to-noise ratio rtOkabe et alJl20l3l . 
We therefore calibrate the isotropic PSF correction, using 
galaxies detected at very high significance, i.e. a signal-to- 
noise ratio of v > 30. This selection acts to suppress the 
measurement uncertainty of (P g ) a p that is caused by low 
signal-to-noise ratio in the objects used for the isotropic 
PSF correction in other studies. The polarizability tensor 
is first estimated by the scalar correction approximation 
( Pg) a p = \tr[Pg\5 a p. We then compute the median of 
( Pg)ap in r g , with an adaptive grid to assemble as uni¬ 
formly as possible, where r g is the Gaussian smoothing 
radius used in the KSB method. We employ a size condition 
of ru > + a-(r’H') and r g > f* + cr(r*) and a positive 

raw P g . Here, f £ (cr(r£)) and f* g ( cj(r * g )) are the median 
(rms dispersion) of half-light radii and Gaussian smoothing 
radii for the stars used for the anisotropic PSF correction 
described above. Although galaxies and stars are well 
separated in rh, we applied the cut in r g so as to exclude 
negative values of P g that are obtained for very small 
galaxies. We also checked the level of stellar contamination 
that galaxy catalogues, selected based on the size cuts 
described here, might suffer. We found that the level of 
contamination is below 1 per cent, mainly due to the fact 
that the number density of faint galaxies is an increasing 
function of apparent magnitude. The faint galaxies therefore 
out far out-number possible stellar contaminants. 

We interpolate the polarizability tensor for individual 
galaxies with v > 10 as a function of r g and the absolute 
value of the ellipticity, \e\. Here, we used |e| instead of each 
component e a because the isotropic PSF correction is per¬ 
formed by the hal f of the trace , -tr\P g \5 a p Then, in a ma¬ 
jor departure from lOkabe et aD (l2013h . we applied a similar 
interpolation for the signal-to-noise ratio, v. An rms error 
of the ellipticity estimates, u g , is estimated from 50 neigh¬ 
bours in the magnitude-r 9 plane. We also experimented with 
smaller and larger numbers of neighbours and found that our 
results are unchanged. 


2.5 Shape measurement tests 


We use two simulated datasets to test the reliability of our 
faint galaxy shape measurements, broadly f ollowing the ap¬ 
proach introduced by the S TEP programme dHevmans et al.l 
120061 ; iMassev et al. l2007tl with important modifications 
compared to STEP and the recent cluster weak-lensing liter¬ 
ature. These modifications are designed to match our data, 
science goals, and our intention to use the weak-shear signal 
on scales of a few hundred kpc to constrain the shape of the 
matter density profile. The first modification is to test our 
ability to measure reduced shears upto g ~ 0.3, as seen in 
the inner regions of clusters. The second modification is to 
produce simulated fits frames that match the angular size of 
the Suprime-Cam data; this allows us to include sufficient 
galaxies with u > 30 that we can test our approach to the 
isotropic PSF correction. This is critical to validating that 


our shape measurement pipeline delivers accurate shapes at 
faint flux levels. We express the results of the tests outlined 
below following the STEP convention of: 

input input , /o\ 

g<x-g a = m a g a v + c a (3) 

where g a and ( 7 )f put are the measured and input elliptici- 
ties, respectively; m a is the multiplicative bias and c a is a 
residual additive term. 

Note that in both of the tests described below, we apply 
a constant shear to all of the simulated galaxies, and thus 
ignore higher order lens ing effects that are present close t o 
the Einstein radius le.g. lOkura, Umetsu &: Futamasell2007l f . 
Higher order effects are negligible on the scale that we mea¬ 
sure and fit the shear profile of clusters. The distribution of 
Einstein radii for a background redshift of 2 = 2 for clus¬ 
ters from our sample that are known stron g lenses is log¬ 
norm al, peaking at #e(Ts = 2) = 14.5 arcsec (IRichard et al.l 
l20ldf . Rescaling this to the typical redshift of z ~ 0.8 for 
the red background galaxies that we use in this analysis, we 
estimate 8e(zs = 0.8) ~ 12 arcsec. Converting to physical 
projected distances, we therefore estimate an upper limit 
of 30/i^ 1 kpc on the typical Einstein radius at the median 
redshift of our cluster sample. Note that we regard this as 
an upper limit because half of our cluster sample have not 
been identified as strong lenses, and thus likely have smaller 
Einste in radii than the clusters discussed by i Richard et all 
(120101) . For comparison, the innermost radius to which we 
typically fit the shear profile in Section 13.21 is 150h 1 kpc. 
Therefore, our shear analysis begins at clustercentric radii 
of > 5 #e, i.e. on scales where higher order lensing effects are 
negligible. 

The first simulated dataset follows lOkabe et al.l (120131 . 
note that these authors also tested their code upto g ~ 0.3), 
and is based on simulated images, kindly provided by M. 
Oguri, that are gene rated with toy models using the software 
STUFF jBertinll2009l 'l. Each galaxy is characterized by bulge 
and disc components, with Sersic profiles indices of n = 4 
and 1, respectively. Galaxy images are convolved with a PSF 
model based on the Moffat profile E(R) oc (l + (R/a) 2 ) ^, 
with seeing in the range 0.5 ^ FWHM Sj larcsec and the 
M offat profile with p ower slopes 3 < /3 < 12, as described 
in lOguri et al] (|2012T ). A number of fits frames matching 
the Suprime-Cam field of view were produced and analysed 
using the pipeline described in i)2.4l We obtain a shear cali¬ 
bration bias of m a ~ —0.02 and c a ~ 10~ 4 (Fig. [2j. 

We extend lOkabe et akl ’s tests with a second simulated 
datas et, using the SHERA software (iMandelbaum et al.l 
120121) to generate simulated ground-based observations that 
match the properties of our observational data. The galax¬ 
ies images included in these simulations are from the COS- 
MOS Hubble Space Telescope observations, as described by 
IMandelbaum et al.l in detail. We convolve the simulated 
data with a Moffat profile that matches our observational 
data: /3 = 4.6 and FWHM = 0.7". We generated both non- 
rotated and 90-degree rotated images to extract a shear es¬ 
timate from a galaxy pair because the intrinsic ellipticity 
cancels out. The magnitude, size, and signal-to-noise ratio 
(u) distributions of the simulated galaxies match those of our 
Subaru data. We analyse these data using the same pipeline 
as above, obtaining again m a ~ —0.02 and c a ~ 10~ 4 - 
i.e. consistent results from two methods of simulating the 
Subaru data (Fig. [2j. 
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Figure 2. Top - Results of the shear measurement tests, showing the multiplicative bias, m a , for simulated observations (0.7 ,/ seeing) of 
model galaxies (GLAFIC; Left) and realistic galaxies (SHERA; Right). The signal-to-noise ratio for the detection are v > 10, matching 
the sample of galaxies that we select for our analysis. Red diamonds, and blue circles, green triangles and magenta squares denote mi and 
m 2 for GLAFIC and mi and m 2 for SHERA, respectively. The lower subpanels show the deviations from the input values. The dashed 
line shows ^output = 9i n put, arid the solid line shows the result of the test described in the Section 2.5, in both upper and lower subpanels. 
Bottom — The shear measurement bias, m a , as a function of the magnitude and the Gaussian size using the imaging simulations of 
realistic galaxies (SHERA). Green triangles and magenta squares denote mi and m 2 , respectively. The bias is not a strong function of 
galaxy size or magnitude. 


We also checked the magnitude and size dependence of 
the shear calibration using the SHERA-based simulations. 
The low level of shear bias detected above does not show a 
strong trend with size and magnitude, with m a ~ —0.02 and 
c a — 10 -4 (Fig. [2]) down to apparent magnitudes of %' = 26. 
This result is achieved because of the high signal-to-noise 
threshold that we apply to galaxies used for the isotropic 
PSF correction described in 42.41 

In summary, our shape measurement bias is below our 
4 per cent goal and does not depend on the size of galaxies. 
However we note that due to the finite number of galaxies 


used in our tests, in particular in the SHERA-based test 
(due to reliance on the COSMOS dataset), we cannot rule 
out the possibility that our shape measurement biases are 
different from those obtained here. 

2.6 Photometry and redshift estimates 

We will select faint background galaxies based on their 
location in the (V — colour magnitude plane in 

42.71 We therefore ana lyse the data using SExtractor 
dBertin fe Arnoutslll996l ). adopting MAG_AUTO as the total 
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i'-band magnitude of each object. (V — i') colours are mea¬ 
sured in seeing matched frames, within an aperture of 1.5x 
the FWHM of point sources in the poorer resolution of the 
two reduced frames for each cluster. Hereafter for conve¬ 
nience we often denote colour as C = (V — i'). 

Early-type cluster galaxies occupy a narrow well-defined 
relation in the colour-magnitude diagram - the so-called red- 
sequence. We fit a linear model of the form 


{V - i') E/S0 = ai'+ b. (4) 


to galaxies at i ^ 22 in order to define the colour of the red 
sequence as a function of i'-band magnitude for the purpose 
of selecting galaxies relative to the red sequence in H2.7I 

Our overall strategy is to combine location in the colour- 
magnitude plane with redshift estimates and reduced shear 
measurements to identify a low contamination sample of 
background galaxies. The next step is therefore to estimate 
the redshift of each galaxy in the photometric catalogues. 
We base these estimates on the COSMOS UltraVISTA pho¬ 
tometric redshift catalogue (McCracken et al. 2012; Ilbert et 
al. 2013) that is limited at i' < 27.5, and benefits from four 
deep near-infrared filters Y, J, H and Ks . This filter cover¬ 
age enables more robust photometric redshifts for galaxies 
at z > 1.3 than earlier versions of the COSMOS catalogue, 
since the Balmer break is redshifted to the near-infrared 
for these galaxies. Furthermore, the COSMOS UltraVISTA 
photometric redshifts are tested against almost 35,000 new 
spectra with galaxies at z > 1.5 (for more details see Ilbert et 
al. 2013). This catalogue provides currently the most reliable 
redshift distribution for a magnitude-limited galaxy sam¬ 
ple that reaches i' ~ 26 with Suprime-Cam on Subaru. We 
emphasise that the COSMOS dataset includes observations 
through the same filters with the same camera mounted on 
the same telescope that we use in this study, thus matching 
LoCuSS and COSMOS photometry is straightforward. 

The lensing kernel for the i-tli galaxy in our photometric 
catalogues, hi = Di s p/D s p, is estimated by an ensemble 
average of the N nearest neighbours in colour-magnitude 
space of the i-th galaxy in the COSMOS catalogue: 

1 N 

hi = (Dis/D s ) cosmos = jy ’^2 Dis,j(zs)/D a j(z s ). (5) 

3 

Here, D s and Di s are the angular diameter distances from 
the observer to the sources and from the lens to the sources, 
respectively. The source redshift, z s , is the median of the 
likelihood distribution for the photometric redshift of each 
COSMOS galaxy. We estimate the uncertainty on hi as the 
sum of the individual photometric errors and standard errors 
of the sample: 


2 


1 N 

= -T 

N ^ 


3* l 

OjgCOSMOSj + 


1 

N - 1 


^COSMOS 



The typical uncertainties are erg ~ 13 — 28%. We adopt N = 
100 , and check that our results are insensitive to whether 
we adopt N = 50 or IV = 200 finding that the redshift 
estimates for individual galaxies are randomly changed by 
a few per cent. We include this small uncertainty in the 
redshift uncertainties that are incorporated in to the error 
bars on the cluster shear profiles in Section m 


2.7 Selection of background galaxies 

Contamination of background galaxy catalogues by unlensed 
member/foreground galaxies leads to a systematic underes¬ 
timation of reduced shear signal. This is often referred to as 
a dilution effect, because the contaminants dilute the signal 
( [Broadhurst et all 120051 : lUmetsu et al.l l20ld : lOkabe et al.1 
l2010al . 20131 '). The dominant source of contaminant galax¬ 
ies is faint cluster members, the number density of which 
increas es towards the clus ter core. 

In lOkabe et al.l (120131 1 we quantified the contamination 
level as a function of the colour offset of faint galaxies from 
the red-sequence, defined by AC = (V— i')— (V— i^bso, and 
concentrated on AC > 0, i.e. galaxies redder than the red 
sequence. We here briefly summarize the method. The mean 
tangential distortion strength is averaged over all galax¬ 
ies satisfying a given colour cut, across all 50 clusters and 
all cluster-centric radii. The mean lensing signal increases 
strongly as a function of AC close to the red sequence, flat¬ 
tening to a shallower trend at larger colour offsets (see Fig. 1 
of Okabe et al. 2013). We interpreted this behaviour with a 
two component model comprising contamination by cluster 
members and the redshift dependence of shear signal - i.e. 
redder galaxies are on average more distant and thus present 
a stronger lensing signal. The model includes a parameter 
that describes the fraction of the total population of galax¬ 
ies that are contaminants from the cluster. Thus by fitting 
this model we were able to measure the contamination level 
and adopt a colour cut that gives a desired level of contami¬ 
nation without relying on any assumptions about the distri¬ 
bution of mass in the clusters, and about the run of number 
density of background galaxies with cluster-centric radius. 
We adopted 1% contamination limit which is less than the 
statistical error of the average mass measurement for the 
sample of 50 clusters. However, this conservative approach 
yields a mean number density of background galaxies of just 
(ribkg) ~ 5.3 ± 1.9arcmin~ 2 , which makes it difficult to mea¬ 
sure individual cluster masses without invoking assumptions 
about the shape of the cluster mass density profile. 

We develop a new method to select securely a larger 
sampl e of background galaxies than achieved bv lOkabe et al.l 
J2013I1 . In the new method, we average the mean tangential 
distortion strength over all galaxies satisfying each colour 
cut (> AC) and all 50 clusters, in several radial bins. The 
stacked lensing strength in the i-th radial bin is estimated 


(A E+>(r0 = 


Y n 5 f +,n5jcr, n^n 


(7) 


where and AE+ is the dimensional reduced shear and the 
subscript n denotes the n-th galaxy located in the an¬ 
nulus spanning n < ri < r 2 . Note that the reduced 
shear AE+ is different fr om the definition of, for example, 
iMandelbaum et al.l d2006l 'l. because we define AE+ in terms 
of the reduced shear, g , and not in terms of shear, 7 . The 
reduced tangential distortion component, g+, with respect 
to the cluster center is defined by 


g+ = ~{gi cos 2 ip + g 2 sin 2 ip ), 


(8) 


where <p is the position angle between the first coordinate 
axis on the sky and the vector connecting the cluster center 
and the galaxy position. The reduced tangential shear is 
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Figure 3. Upper left: The mean lensing signal for all 50 clusters as a function of the colour offset AC in cluster-centric radial bins. The 
innermost bin has the largest shear signal (upper most set of points overlaid with brown curve), with progressively smaller shear signals 
seen at larger cluster-centric radii (lower sets of points overlaid with red through yellow, green and blue curves). All of the curves are flat 
at AC > 0.5 - i.e. at galaxy colours used for our analysis in Okabe et al. (2013). At smaller values of AC the curves relating to smaller 
cluster-centric radii steepen whilst those relating to larger cluster-centric radii are much flatter. The best-fit models are represented by 
the solid coloured curves. The legend states the typical weighted harmonic mean clustercentric radius of each radial bin. The dashed 
magenta curve denotes the location of the 1% dilution level; the intersection between this curve and the solid curves indicates the colour 
cut required in each radial bin to achieve 1% dilution. The vertical dotted green line is the lcr width of the red sequence determined 
by the bright member galaxies. Lower left: The best-fit models of the dilution level of member galaxies as a function of the colour 
offset. The colours are the same as those in the upper left panel. The dilution level at innermost bin (brown) shows ~ 40% at AC = 0. 
The horizontal dashed line (magenta) is our requirement for the dilution level. Right: The colour-cut for background selection with 1% 
dilution as a function of clustercentric radius. The colour cut is AC ~ 0.5 at r ~ 100 h~ 1 kpc and AC = 0 beyond ^ 1.3 h~ 1 Mpc, 
respectively. 


expressed by 
7+ 


9+ = 


1 — K, ’ 


(9) 


where 7 + is the tangential shear and k is the dimensionless 
surface mass density. The weighting functions are taken into 
account for both the statistical weights, w„, and the criti¬ 
cal projected mass density, S cr . n, describ in g the lensing effi¬ 
ciency. The weighting w n (e.g. Hoekstra^ftam^&^J^uijken 


et al.l2Q0l:lOkabe et al.l2010al:IUmetsu et al 


200d:lH 


20ld: lOeuri et al. l2012h is used to down-weight galaxies 
whose shapes are less reliably measured, based on the un¬ 
certainty in the shape measurement, 07 , n , given by 


ry2 1 q -2 Y? 

1 ^ g,n ^-'cr,n 


( 10 ) 


We choose a = 0.4 throughout this paper. The critical pro¬ 
jected mass density for individual clusters is defined by 


Yj cr — 


47 tGDi 


( 11 ) 


where D; is the angular diameter distance to a cluster and 
f) is the lensing depth (equation 0 . 

To clarify the relationship between AE + and g + in the 
presence of contamination of background galaxy catalogues, 
we write an equation to describe the contamination effect, 
as follows: 


(AE+) 


cr,7iU7bkg,n T 0 X W no n —bkg,n) 

X0(u7bkg,Ti T U.’non — bkg.ri ) 

V) AE-|_,nU7bkg,n 
/(D(^bkg ,n T ^non-bkg,n) 


( 12 ) 


where AE+ = g+ E cr is the dimensional shear for each back¬ 
ground galaxy and the subscripts “bkg” and“non-bkg” de¬ 
note background galaxies and non-background galaxies, re¬ 
spectively. Since non-background galaxies are not lensed by 
the cluster, the second term in the numerator is zero. There¬ 
fore, contamination by non-background galaxies causes the 
observed lensing signal to be underestimated. The underesti¬ 
mation of the ensemble average is simply determined by the 
fraction of non-background galaxies (mainly member galax- 
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10 Okabe and Smith 


ies) to background galaxies. This is even clearer if, for the 
sake of illustration, one assumes that the galaxies have uni¬ 
form weights, w = 1, in which case Equation [12] can be 
rewritten as: 

1 (13) 


(AS+) = (AE + )bk g 


1 T n non —bkg/^bkg 


where (AE+)bk g is the average distortion strength using a 
pure sample of background galaxies. The correction factor, 
(1 + n non _bkg/nbkg) is equivalent to the boost factor, under 
the assumption of a radially uniform distr i bution of back¬ 
groun d galaxies (e.g. lAppleeate et al.ll2014l ; iHoekstra et al.l 
120151). that we discuss in Section [4.1.41 Equation 1121 shows 
that the average lensing signal obtained from our formalism 
is simply underestimated by the fraction of non-background 
to background galaxies. Therefore also note that the non¬ 
linear term in the reduced shear g+ = 7 +/(l — «) oc AE+ 
does not impact on our estimated levels of contamination. 

The stacked lensing signal is a decreasing function 
of clustercentric radius and an increasing f unction of AC 
for small clustercentric radii, reminiscent of lOkabe et al.l ’s 
(2013) analysis of all clustercentric radii as a single bin (up¬ 
per left panel of Fig. n. However the increase in lensing 
signal at moderate values of AC becomes progressively less 
pronounced as one considers radial bins at larger cluster¬ 
centric radii. This is qualitatively consistent with faint clus¬ 
ter galaxies being the dominant source of contamination, 
given that the number density of galaxies in clusters is a 
declining function of clustercentric radius. 

To interpret quantitatively the mean lensing strength 
as a function of cluster-centric radius, we parameterise the 
galaxy distribution in terms of the projected clustercentric 
radius r and the colour offset AC, as follows: 


n(r, AC) = n b (r) [1 + Bf(AC)n m (r )]. 


(14) 


where n b (r) is the radial distribution of background galax¬ 
ies. We do not assume any specific functions of n b be¬ 
cause the background distribution is not constant, and 
may be depleted or boosted by a magnification bias (e.p 


Broadhurst. Taylor fe Peacock! fl995l: lUmetsu et alj 1201 lL 
20141 ; Coupon. Broadhurst fe Umetsul 2013lb The second 


term in the bracket denotes the member galaxy distribution; 
B is the fraction, /(AC) is the colour distribution and n m (r) 
is the radial distribution. The effective lensing strengths for 
red galaxies are obtained by integrating over the projected 
radius and the colour offset: 


intrinsic scatter in the colour distribution and the photomet¬ 
ric error. We assume that er is radius-independent. When we 
av erage equatio n (11511 over all radial bins, the formulation 
in lOkabe et al.l ( 20131 ) is recovered: 


«AE+» = ^(AE+,4 ~ A [l - BF(AC)] . 


(17) 


(AE + ,i)(n <n < r 2 , AC <) = 


f? c d(AC)f, 


r 

JL 


A C d ( AC )l 
At 


1 + BF(AC)Ni ' 

Here, Ai = AE +09 is the lensing signal estimated from the 
pure background galaxies and is thereby determined by the 
cluster mass distribution. The contamination levels in the 
colour and radial distributions are described by F(AC) = 
d(AC)f(A C) and Ni = f’/ 1 drn m (r)r, respectively. As 
in lOkabe et al.l l|2013h . we employ a Gaussian distribution 
centering at AC = 0 as the colour distribution of member 
galaxies, 

F(AC) = [1 - erf(AC/%/2cr)]/2, (16) 

where a is the width of colour distribution composed of the 


In this paper, we simultaneously take into account the colour 
distribution and the radial distribution for member galaxies. 
We employ n m (r) = exp ( — r/rp) as the r adial distribution of 
member galaxies ((Applegate et al]|2014f ). Subsequently, fit¬ 
ting is performed to obtain the colour and radial distribution 
of member galaxies. 

In summary, the fitting parameters are Ai, B, a and 
r 0 . We stress that this method does not assume any specific 
mass models, which is important to interpret the results af¬ 
ter defining the background sample. As shown in the upper 
left panel of Figure^] the best-fit model (solid lines with dif¬ 
ferent colours) well describes the data. The best-fit colour 
width, a = 0 . 21 , is higher than the mean width expected 
from the intrinsic scatter determ ined by the cl uster bright 
galaxies, which is consistent with lOkabe et alj (l2013l l . This 
large value of a is driven by the statistical scatter of the faint 
galaxies included in the calculation - i.e. the photometric 
uncertainties at i' ~ 25. We note that our method assumes 
that the colour distribution of galaxies redward of the red 
sequence is Gaussian; the large value of a therefore helps to 
ammeliorate any concerns that the wings of the actual distri¬ 
bution contain an excess of galaxies over the assumed Gaus¬ 
sian form. The characteristic radius of the member galaxy 
distribution is r 0 = 258 h _ 1 kpc. As expected based on 
the previous qualitative discussion, the highest level of con¬ 
tamination occurs at AG = 0 at the smallest clustercentric 
radii, with F ~ 0.4 (Fig. [3]), and the level of contamination 
declines significantly with increasing clustercentric radius. 
Contamination is negligible in the cluster outskirts. We con¬ 
servatively adopt a limit of 1 % on contaminating fraction, 
and use this to define a radially dependent colour cut (right 
panel of Fig. [3]). Note that at r > 1.3 h _1 Mpc the con¬ 
tamination level is so low that we adopt AC > 0 in this 
region. We achieve a number density of background galax¬ 
ies of ribkg — 5 — 20 arcmin ~ 2 (Table [TJ, with a mean of 
(n-bkg) — 12 . 8 arcm in ~ 2 that is more than double that of 
lOkabe et al.l (l2013l l. 

We also use our mid- and far-infrared observations with 
Spitzer and Herschel as a sanity check on the possible impact 
drAFj+n b {r)r 0 f heavily dust-obscured galaxies on our red galaxy selec- 
|r tion and on our COSMOS-based estimates of /3 in the previ¬ 
ous Section. Specifically, we consider whether dusty cluster 
(15)members might leak into the red background galaxy sam¬ 
ples and whether /3 might be biased due to the presence 
of optically faint and heavily dust-obscured galaxies - i.e. 
Luminous Infrared Galaxies (LIRGs) and Ultra-Luminous 
Infrared Galaxies (ULIRGs). On the latter point, the COS¬ 
MOS photometry extends to 8p,m, whereas we have ob¬ 
served half of the cluster sa mple discussed in t his article 
with Spitzer/MIPS at 24/rm (lHaines et a.l .1 EfoTHli and with 
Herschel/PACS and SPIRE at 100 — 500^im ( Smith et alJ 
l 20 ldh albeit only to a depth corresponding to a bolometric 
infrared luminosity of Lir ~ 5 x lO lo Z /0 at 2 ~ 0.2. The 
spectroscop ic completeness of follow-up observations with 
Hectospec dFabricant et al.ll2005ll is 96 per cent for objects 


drn(r 
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Weak-lensing mass calibration of galaxy clusters 11 


detected with Spitz er down to 0.5mJy and and ~ 80 pe r 
cent down to i — 20 dMulrov et al.ll2014l : lHaines et al.ll2015lj . 
Clearly these data are not sensitive enough to obtain defini¬ 
tive estimates of the number and redshifts of dusty galaxies 
that satisfy our faint red optical background galaxy selec¬ 
tion. However they provide a useful sanity check based on 
bright gal axies. We select gala xies from the catalogues dis¬ 
cussed bv lHaines et al.l J2015l l at i < 20 that have {V — i) 
colours that would place them in our red background galaxy 
catalogues if they were faint enough. We find that ~ 1 per 
cent of these bright optically red galaxies are LIRGs, and 
~ 1 per cent of the same bright optically red galaxies are 
cluster members, and there is no overlap between these two 
populations. We therefore conclude that the bright optically 
red galaxy population seen along lines of sight through our 
cluster sample appear to be consistent with (1) LIRGs and 
ULIRGs not being a significant population in our red back¬ 
ground galaxy samples, and thus not being a concern in 
terms of the accuracy of /3, and (2) our red background 
galaxy catalogues suffering just 1 per cent contamination 
by cluster members. 


3 MODELLING AND RESULTS 
3.1 Model Fitting Methods 


mass model, 


/model(ri) 


AE mo del {Ti ) 

1 R/E mo dei(i7)) 


(19) 


where E = E ci k, and AE = E cr 7 are the convergence and 
the shear in physical units, respectively, and k and 7 are 
the dimensionless convergence and shear, respectively. The 
factor, Ki for the i-th bin is given by: 


K, = 


E.» El cr,nWn 


( 20 ) 


Note that Ki is computed separately for each bin due to the 
radial dependence of the redshift of the background galaxies. 
The covariance matrix, C, in equation M is given by: 


C = C g + Cs + Clss- 


( 21 ) 


Where the shape noise, C g , in each radial bin is estimated 
as 


Cg,ij — 


1 2 v-i2 .2 

1 / j n £h?,n^-'cr,n^n 

2 (E„«M 2 


( 22 ) 


where 5ij is the Kronecker delta and the factor of 1/2 ac¬ 
counts for the rms noise, <J g: n, of two distortion compo¬ 
nents. The photometric redshift error matrix, C s , is com¬ 
puted from: 


Cs,ij — Cs 


AElmod el Em 0 d el cr ir,i 
(1-Apodal(r;)) 4 


(23) 


We describe how we compute the reduced tangential shear 
profile of each cluster and the model fitting procedure. We 
apply the methods described in this section to measure clus¬ 
ter masses in the next section. 

We centre each cluster shear profile on the centroid 
of the optical emission from the cluster’s BCG, following 
numerous previous studies that have shown BCGs to be 
a reliable cluster centre for w eak-lensing s t udies o f mas ¬ 


sive galaxy clusters, including IOkabe_et~al. (2010a, 


Ivon der Linden et al.l ll20l4l . and Hoekstra et al 


Iml, 


20151) . 


For example, we derived an upper limit of 32/U 1 kpc on 
the mean offset of BCGs from the underlying centre of the 
cluster mass distr ibution for the sample studied here, in 
lOkabe et alj (l2013ll . This upper limit is a factor of 5 smaller 
than the typical innermost radius of the shear profiles upon 
which our mass measurements are ultimately based in Sec¬ 
tion 13.21 Any bias caused by centring our shear profiles on 
the BCGs is therefore negligible. 

The reduced shear in a given annulus centered on a 
given cluster is computed by azimuthally averaging the mea¬ 
sured galaxy cllipticities, as defined by Equation [3 The 
mean redshift of the background galaxies is a function of 
cluster centric radius, due to our radially-dependent colour 
cut (Section 12.71) . The formulation in equation [3 takes ac¬ 
count of these differences by expressing the reduced shear in 
physics units. 

We employ a maximum-likelihood method to model the 
shear profiles, and write the log-likeliliood as follows: 


— 2 In C = ln(det(CV,)) + 


where the first term is given by 

+ _ 2AE-|_) 2 W)^((Ti; i „/Ec ri n) 2 J 

Ls.tJ — 7T=; 777 Oij. 

V (E„»J 2 J 

The second term of the equation (1231) is the photometric red¬ 
shift errors through an error of the conversion factor, gk, 
in the mass model ©• The covariance matrix of uncor- 
relate d large-scale structure (LSS), Clss, along the line-of- 
sight llSchneider et al.l 1 199a) at an angular separation be¬ 
tween 9i = n/Di and 0j is given by 

C^ss.ij = j l ^P.(l)J2{Wi)J 2 (Wj), (24) 

where P K (l) is the weak-lensing p ower spectrum (e.g. 
ISchneider et al.1 1 19981 : IHoekstral 120031 ') , calculated by mul- 
tipole l, the source redshift, and a given cosmology. We 
empl o y the redshift, m in (z s ,j, z s ,-j), and WMAP9 cosmology 
(jHinshaw et al.ll2013h . And, J 2 (l9i) is the Bessel f unction of 
the fi rst kind and second order at the i-th annulus (IHoekstral 
l2003f) . 

It is also important to compute the radius of each radial 
bin correctly, because systematic errors in the placement of 
the binned shear measurements on the radial axis can cause 
systematic errors in the mass measurement when a model 
is fitted. This is particularly important in practice, because 
the number density of background galaxies is neither uni¬ 
form nor infinite. As described in detail in Appendix m we 
found that the best radius at which to place the measure¬ 
ment of mean tangential shear in a radial bin is the weighted 
(18) harmonic mean: 


^ /model { r i ) ) Cij (AE-j-J f model(u)) 5 /j.\. E n Wn 

hi 1 E„ 1 ’ 


(25) 


where the subscripts i and j are the i— and j —th radial bins. where w n is given by Eauation llOl We therefore compute bin 

Here, /model is the reduced shear prediction for a specific radii in this way, in physics units, in the rest of our analysis. 
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12 Okabe and Smith 


Finally, before fitting models to the shear profiles (Sec¬ 
tion [3j2]), we calculate the signal-to-noise ratio of the tan¬ 
gential shear profile as follows: 

(. S/N f = Y, AS+,<(£?,,y + C sAj + C Lss , ij )“ 1 AE +J , (26) 

ij 

Each cluster is detected individually at signal-to-noise ratio 
of 3 < S/N < 11 (Table 0. 


3.2 Mass Measurements 


To infer galaxy cluster masses from the shear profiles, we fit a 
model to the latter. For this purpose we adopt the universal 
mass density profile dNavarro, Frenk fc Whitdll996l . Il997bl . 
hereafter NFW), that has had considerable success in de¬ 
scribing dark matter halo profile spanning a wide mass range 
by numerical simulations based on the CDM model of struc¬ 
ture formation. It has also been shown that the ensemble 
mass of a sample of clusters can be recovered to good preci¬ 
sion from this approach, provided sufficient care is taken over 
the radial range over which this model is fitted to data ( e.g. 
iBecker fc Kravtsovll201l1 : iBahe. McCarthy fe Kinsll2012i . 

The NFW profile is expressed in the form: 


pNFw(r) 


_ Pa _ 

(r/r s )(l + r/r s ) 2 ’ 


(27) 


where p a is the central density parameter and r s is the scale 
radius. The three-dimensional spherical mass, Ma, enclosed 
by the radius, r/\, inside of which the mean density is A 
times the critical mass density, p cr (z), at the redshift, z, is 
given by 


Mn 


r(< r &) = 


47T p s r\ 

c 3 

C A 


ln(l + ca) — 


CA 


(28) 


1 + ca _ 

The NFW profile is fully specified by two parameters: Ma 
and the halo concentration ca = ca/tv We fit this model 
to the shear profile of each cluster, taking full account of 
errors of shape measurements, photometric redshifts and the 
uncorrelated LSS fSection 13.11) . For a given Ma and ca we 
predict the observed shear sign a l follo wing the formalism 
described bv I Wright, fc Brainerdl l|200flh . 

Measurements of Ma are mainly s ensitive to the lensin g 
signal around the overdensity radii r a (lOkabe et al]|2010al l . 
However the concentration parameter, ca, is more strongly 
affected by the lensing signal in the cluster central regions, 
i.e. cluster centric radii of hundreds of kpc. Our careful se¬ 
lection of background galaxies ensures that contamination of 
our background galaxy samples is negligible across the full 
radial range of our shear profiles. However the very stringent 
colour cut employed in the central regions, A C > 0.52, and 
the relatively small solid angle subtended by these innermost 
bins, render them the noisiest of the entire radial range. To 
guard against obtaining results on concentration that suffer 
biases due to the noisy inner profiles, we choose a binning 
scheme for each cluster via the following procedure. We fit 
the NFW model to a suite of measured shear profiles that 
span inner radii in the range n n = 50 — 300 /i _1 kpc, outer 
radii in the range r ou t = 2000 — 3000 h, -1 kpc, and number of 
bins in the range iVbm =4 — 8. We then compute the mean 
of the suite of Ma values obtained from these fits, and adopt 
the binning scheme that yields the value of Ma closest to 
that mean. Note that we allow the virial concentration pa¬ 
rameter to be in the range 0 < c v ir < 30 in the fits. Also, we 


restrict the radial range of the shear profile fits for A1758N 
to r ou t < 2100/i _1 kpc to avoid contamination of lensing sig¬ 
nal by its neighbour A1758S. Note that we test the proce¬ 
dure described above using mock observations of simulated 
and toy model clusters, and confirm that it returns masses 
and concentrations with negligible bias (Section 14.1.11) . 

Figure 0 shows the tangential distortion profiles as a 
function of the projected cluster-centric radius for two ex¬ 
ample clusters, ABELL 2390 and ABELL 0901. The former 
is among the most massive in the sample and the latter 
among the least massive. The tangential shear clearly de¬ 
creases from the cluster centre to the outskirts, with the less 
massive cluster, ABELL 0901, presenting an overall shear 
signal of approximately half that of the more massive clus¬ 
ter, ABELL 2390. Note that the 45 degree rotated compo¬ 
nent times the clustercentric radius r, rAE x , is consistent 
with zero - i.e. this simple test of residual systematics is 
consistent with zero. 

Tabic [2] lists Ma from our weak-lensing analysis, de¬ 
fined as Ma = A p cr i t (,z)47rr^/3 where p cr it ( z ) is the 
critical density of the universe at the respective cluster red- 
shifts, and A = A v ir, 200, 500, 1000, and 2500. We also list 
Ma defined as Ma = Ap m (z)47rr^/3, where p m (z) =p C rit 
(z)fiM(z) is the mean matter density of the universe, and 
A = 180 and 200. We denote these latter two masses as 
Misom and M^oom respectively. 


3.3 Mass-Concentration Relation 


Numerical simulations (e. 

s. Bullock et al. 2001: Duffv et al.l 

20081: Bhattacharva et al. 

20131: iDiemer & Kravtsov 20141: 

Meneehetti et al. 2014!: Ludlow et al.| 201411 predict that the 


halo concentration ca and the mass Ma for the NFW mass 
model is weakly anti-correlated. Such a correlation is natu¬ 
rally explained by the hierarchical structure formation, that 
is, less massive halos first form and more massive halos form 
through mass accretion and mergers of smaller objects. The 
characteristic central density of more massive halos is lower 
as reflected by the critical mass density of the universe at the 
redshift of collapse. Measurements of cluster mass and con¬ 
centration therefore provide us with a unique opportunity 
to test structure formation. 

Our cluster sample is selected purely on X-ray lumi¬ 
nosities without imposing any requirement on the physical 
properties of the clusters. In particular, we do not select on 
the dynamical state of clusters as inferred from their X-ray 
morphology. We are therefore able to investigate the corre¬ 
lation between mass and concentration for a large sample 
of clusters that is unbiased beyond that which is inherent 
to an X-ray selection. A typical cluster in our sample has 
a concentration of C200 — 4 (Figure [5]), with central values 
of C200 in the range C200 ~ 2 — 20. We quantify the mass- 
concentration correlation with the following function: 

c.(M.), m ( m , 4 ^ 1 m J , (») 

where Co and b are the normalization of the concentration 
parameter at Ma = 10 14 h -1 Mg and the slope, respectively. 
This form is motivated by the stu dies of the numerical sim¬ 
ulations fe.g. iBullock et ahll200ll l. Note that we ignore red¬ 
shift evolution in this model because the redshift range of 
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r [h 1 Mpc] r [h 1 Mpc] 


Figure 4. Radial profiles of the tangential shear component (top panel), AE_|_, and the product of the 45 degree rotated component, 
AE X , and clustercentric radius, r 

(bottom panel), , for ABELL 2390 (left) and ABELL 0901 (right), respectively. 


our sample is narrow. When we fit this model, we take ac¬ 
count of the correlation between the errors on concentration 
and mass by calculating the error covariance matrix, and in¬ 
clude the intrinsic scatter of the lensing-based concentration 
parameter, <Ji n t. The log-likelihood is given by: 


— 21n£ 


( a ?nc,i + b 2 a[ 

n M,i 2bai n 

c,ln M,i + *L) 
i 


+E 


(ln(ci) — (a + b ln(Mi))) 2 
°fn c,i + b2(J ?n M,i - 26(71nc,In M,i + + nt ’ 


where a = log(co) and <ri nc and m n M are the fractional er¬ 
rors of the concentration and the mass, (Tinc.inM is the error 
correlation, and (Tint is the intrinsic scatter in Inc. We per¬ 
form the fitting at overdensities of A = 180m, 200m, vir, and 
200. The normalisation of our best-fit mass-concentration 
relation is in excellent agreement with the results of re¬ 
cent numerical simulations at zi = 0.23 (Table [3] Fig¬ 
ure [5] Bhuttadmr va^et al.l 120131 ; iDiemer fc KravtsoyTteOl^ : 
iMeneghetti et al.l l2014l) .~~A.s an aside, we note that these 
three recent independent theoretical studies agree both with 
each other and with our observational results, whilst older 
simulations showed considerable variation between their re¬ 
spective mass-conc entration relations and a lower overa ll 
normalization (e.g. iDuffv et al.l 120081 ; IStanek et al.l I2OI0I) . 
The best-fit slopes agree with the weak-mass dependence of 
the concentration, b ~ —0.1, seen in simulations (Table [3]), 
although the uncertainties are too large to rule out positive 
values of b. Adding less massive clusters and increasing the 
number density of background galaxies will allow improved 
constraints in future studies. 


3.4 Stacked Lensing Analysis 

Stacked lensing analysis is a powerful technique for mea¬ 
suring the average density profile of a sample of clus¬ 


Table 3. Best-fit parameters for the mass concentration relation. 


A 

Co 

b 

<7int 

180m 

Q 04+ 3 " 

y,Z4 -2.50 

-0.20i°+ 

< 0.17 

200m 

0 yyj+3.75 

0. <4_ 2 49 

-0.19±°i| 

< 0.17 

Virial 

* - zd _ 2.07 

—0 17+ 015 
U * A '-0.16 

< 0.18 

200 

c 1 0+2.08 
—1.44 

—0 14"*" 016 
Ui 14 -0.16 

< 0.20 


ters. Stacking the shear signal from a sample of clus¬ 
ters averages over the distribution of internal structures 
and halo triaxiality, and thus overcomes the structural bi- 


ments (e.g. 

Mandelbaum et al. 20061: Johnston et al.l 20071: 

Okabe et al. 

2010a, 

120131: llJmetsu et al.l 201ll. 20141. 2015bl: 

Oauri et al. 

20121: 

Leauthaud et al.l 20121: iMivatake et al.l 

20131: Niikura et al. 

ISnsI). 


We compute the average lensing signal in physical 
length unit centered on the respective BCGs. Note that our 
redshift range is narrow, and therefore the results described 
below are unchanged if we instead use comoving length 
units. Moreover, we have previously tested that adopting 
physical length units, and not scaling length to an overden¬ 
sity radius, yields an unbias ed measurement o f the stacked 
shear profile of our sample dOkabe et al.|[2013h . The inner¬ 
most radius of the stacked shear profile is that at which the 
innermost bin of the stacked profile contains a minimum of 
one background galaxy from each cluster. The outermost 
radius of the stacked profile is the median of the maximum 
physical scale on which the field of view of the Subaru ob¬ 
servations fully encloses a circular aperture centered on each 
BCG. Note that this simultaneously matches the angular 
extent of the data, and satisfies the requirement placed on 
the innermost radius. The stacked shear profile decreases 
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14 Okabe and Smith 


Table 2. Mass measures for individual clusters on (M^) cog 


Name 

M vir 

M 200 

M 500 

Afiooo 

M 2500 

M 180 m 

Af200m 


10 14 /l _1 Mg 

10 14 h- 1 M e 

10 14 h~ x MQ 

10 14 h- 1 M Q 

10 14 h~ 1 MQ 

to 14 h- 1 M Q 

10 14 h~ 1 M(r ) 

ABELL2697 

10.37±:?;?? 

8.0311;“ 

4 62 +0 ' 84 

^• DZ —0.84 

9 79 +O .73 
z - <z -0.80 

1 19 +O. 6 I 

X - 1 Z —0.53 

12-291?;?? 

11.761?;?? 

ABELL0068 

7.78 ±\™ 

6.65li;35 

4 78+ 0 - 78 
<o_0.71 

q ci +0.50 
°' o± -0.48 

9 1 9 + 0.34 
z ' xz —0.36 

8.641??? 

8 . 39 !??° 

ABELL2813 

9 47 + 2 - 53 
- 2.02 

8.17li:® 4 

5 90"*” 1,06 

°- yu —0.97 

4 37 + 0 - 72 
—0.70 

2 66”*”°' 58 
z.oO-o gs 

10.431?;?? 

10 14+ 2-87 
1u ' 4 ^-2.25 

ABELL0115 

9 49 + 4,56 
.94 

7 04+2-66 
' * u ^-1.97 

q 77+1-14 

2 07+0-77 
z,u -0.77 

0 75+0-51 

U ‘' °—0.40 

ii-eol?;?? 

11-06!?;?? 

ABELL0141 

5.67t 4 ; 4 ® 

4 71+108 

4 —0.95 

3 49 + 0.65 
xy —0.60 

9 99 +O .48 
^•^^_0.48 

1 91 +0.38 
x,zx -0.38 

6-421?;?? 

6 - 22 !?;?? 

ZwC10104.4+0048 

9 7-1 + 1.66 

z -' 4 —1.09 

9 24+1-01 

1-471®;?? 

q 99+0 48 

u.y »_ 0 70 

0 51+°- 46 

0.53 

3.081?;?? 

2.98!?;?? 

ABELL0209 

15 42+ 3 14 

12-75!?;?? 

8-641?;?? 

5 99 +O 68 

3 26”*” 0,46 

o.zu-o 49 

1 7 C7+3.90 

4 ‘- 0 '-3.07 

17.011?;™ 

ABELL0267 

7 90 +I .58 
' *^0-1.33 

5-97li;Sl 

3 92+ 0,64 

o.»^_0.59 

2 63+ 0,43 

Z.OO _ 0 42 

1 q5+0-30 
I.OO-Q 31 

Q qq+1.95 

o.oo-i go 

8.04!!;?? 

ABELL0291 

7 47 + 2.68 

—1.92 

5-63141 s 

q 1 9+0.71 
O,xz -0.67 

1 77 +O. 5 O 
' ' -0.52 

0 - 68 lg-Jg 

9.051?;?? 

8.64!?;?? 

ABELL0383 

6 27"*" 1,83 
°- z '-1.42 

c 9q+l-30 
°- zo -1.07 

3 54+0-65 
-0.60 

2 59 +O 39 

^.O»_o 37 

i-48lg;S 

7 11 +2.30 

1 • i± -1.71 

6.90!?;?? 

ABELL0521 

6 76 +1,58 
°- '°-1.36 

5-6i!i;Sl 

0 77 +O .69 

01 1 -0.65 

2 59 + 0 - 51 
z.oy _ 0 51 

1 39 + 0 - 40 
x .o»_o 40 

7.651?;?? 

7 40"*" 1,83 

'•^ u -1.54 

ABELL0586 

7-69t^ 9 5 

6-65l?;Jl 

5 04"*" 1,12 

0 ,u —0.98 

q Q9+0-69 
o-»^_o.64 

2 60”*”°' 53 
Z.OU-o 68 

8.531?;?? 

0 39+3.54 
°' oz —2.32 

ABELL0611 

10.891??? 

9 23”*” 1 ' 92 
y * zo -1.71 

5 37+1-17 

D *°'-1.09 

4 54 + 0 - 82 
^• O1 -0.80 

2 54 +°- 58 
^•^-0.60 

12 14+ 2-91 
4 Z - 4 ^_2.45 

11-771?;?? 

ABELL0697 

12.62!?;?? 

9-74l 2 2 ;?° 

5 39+ 1 - 08 
0.017—1 00 

3-051°;?? 

1 1 7 +O .56 

4-4 ‘ —0.51 

14 91+ 6 - 72 

14-22!°;?? 

ZwC10857.9+2107 

0 1 9+1.58 
°‘ 4 Z — 1.20 

2.45+°? 

1 45+0-69 
J-. 4 O -0 75 

0 - 88 l®;°? 

0 37+0 - 45 
u -° -0.42 

3.671??? 

q 59 + 1-97 

O.OZ_i 39 

ABELL0750 

7 66”*” 3,98 

( 00_ 2 34 

6.30+71 

4 04 + 1.20 
^• ox —0.95 

q r» 9 + 0.60 

O.U Z _ 0 54 

1 67+0- 39 

4 '°' —0.46 

0 70 + 5 .I 4 
*'°-2.85 

8 . 51 !?;®° 

ABELL0773 

n - 11 till 

9.56+?® 

7 05”*" 0,75 

' .UO_Q >70 

c qq+0.50 
O.OO-Q 48 

q q 7 + 0.37 

0,0 -0.39 

1 9 qi + 2.00 
iZ.Oi-i 70 

12 . 00 !!;®! 

ABELL0781 

8-17 ±1;“| 

6 57 + 1 97 

°*°'-1.65 

3 97 + 1,11 
-1.19 

2 46+ 0,93 
z,ftD - 1.10 

1 09”*"°' 77 
x,uy —0.69 

9.40l?;®° 

9 02 ”*" 3,72 
y,uz -2.60 

ZwC10949.6+5207 

5-o2tlS 

4 44+ 132 
^•^- 1.10 

3 48+0- 79 
0'^°—0.73 

9 7Q+0-56 

z -' 0.53 

1 96”*” 0,44 

4 .yo_o 53 

5.47l?;®° 

5 35 + 1.97 

O.oO-i 49 

ABELL0901 

0 11+1-26 
°' xx —0.97 

2.651+? 

1 95 + 0 - 57 
i.yo_o 50 

1 48”*” 0,39 

.3 6 

0 93 + 0,29 
u.»o_o 31 

3.481?-°? 

3 39+ 147 

o.o»_i 10 

ABELL0907 

18.97±|;“ 

14.281+® 

8.071??? 
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smoothly as a function of clustercentric radius (Figure [U]), 
and yields a signal-to-noise ratio of is S/N ~ 35.6, after 
taking into account the LSS covariance matrix, Clss- 

To interpret the average mass profile from the stacked 


lensing signals, we consider three mass components, 

ASmodel = AEpt + AEnFW + AE2h, (30) 

where AE p t is a point mass associated with the BCGs, 
AEnfw is the large-scale cluster mass distribution that 
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Figure 5. The observed distribution of the concentration param¬ 
eters C 200 as a function of the cluster masses M 200 for 50 clusters. 
The errors denote 68% confidence intervals. The thick and thin 
lines (red) are the best-fit function and the errors, respectively. 
The dashed blue, dotted green and dotted-dashed magenta lines 
are the mean m ass-concentration re l ation from recent numeri¬ 
cal sim ulat ions o f Bhattacharva e t al.1 (l2013l) , rDiemer &; Kravtsovl 
(|2014h and iMeneghetti et al nm at zi = 0.23, respectively. 


Figure 7. Confidence intervals on mass and concentration from 
stacked lensing analysis. From left to right, the contours show 
the low-mass sample (thin green), the full sample (thick red), 
and the high-mass sample (thin blue). Crosses denote the best-fit 
parameters and the contours show the 68.3%, 95.4%, and 99.7% 
confidence levels. The dashed, dotted and dashed-dotted lines are 
the same as Figure [5] 



0.1 1 
r [h~ l Mpc] 


Figure 6. Stacked tangential shear profile for 50 clusters. The er¬ 
rors are composed of (C g ,ti + C s ,a + Clss u) 1 ^ 2 • Thick solid red 
and dashed green lines are the total mass and the NFW model, re¬ 
spectively. The dotted blue and dashed-dotted magenta lines are 
the two-halo term and point source multiplied by 10, respectively. 


we parametetr ise following NFW, and A£2h is the two- 
halo term (e.g. Johnston et al. I l2007l: lOeuri fc Takadall201ll : 


lOgiiri fc Hamanall201ll) to account for structure adjacent to 

the clusters. Note that the latter two terms were ignored in 
the modeling of individual clusters because the noise level 
in individual cluster shear profiles renders them insensitive 
to these contributions. 

We describe the contribution from the point mass, of 
mass M p t, as: 


AEpt — 


Mr 


P t 


(31) 


and adopt a prior on M p t based on the stellar mass for 
the BCG. The stellar mass of each BCG is es timated from 
the A-band luminosity with a ISalpeted (1 19551 1 initial mass 
function. The prior on the point mass then matches the 
mean and standard deviation of the BCG stellar masses. 
The two-halo term is com puted following the formulation of 
lOguri fc Hamanal (l201lli . We use the WMAP9 cosmology 
dHinshaw et al.l 20131 1 to compute the linear power spec¬ 
trum. Given the average mass and redshift for an ensem¬ 
ble of clusters, AE2 h is proportional to pm{zi)bh(M), where 
bh{M) is the halo b ias. To estimat e bh(M), we use a single 
scaling relation dTinker et al.|[2Qlol l which is calibrated by a 
large set of numerical simulations. 

The best-fit model describes the data very well (Fig¬ 
ure El Table [p. The two-halo term is an order of magnitude 
less than the NFW model, with an estimated halo bias of 
bh{M 2 oo) — 5.5. The point mass is constrained by the up¬ 
per limit, Mpt < 6.19 x 10 11 h _1 Mg. The mass and the 
concentration at A = 200 (Figure Q i s in excellent agree- 
ment with both numerical simulations (IBhattacharva et al.l 
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16 Okabe and Smith 


Table 4. Best-fit parameters for stacked lensing analysis for 50 
clusters. 


A 

m a 

CA 


10 14 h~ 1 MQ 


180m 

8-66t°;" 

k e:7+0- 3 7 
°- 0 '-0.35 

200m 

8.39t°;« 

c qo+0.35 
0.33 

vir 

7 65+ 0 - 38 
'• oo -0.36 

4 .D< 5_0 30 

200 

c o't-I-0.28 

U-0< -0.27 

3 69"*” 0 ' 26 

O.u»_o 24 


straints on M 200 and C200 are in excellent agreement with the 
NFW-based measurements (Table [4j, and also agree within 
~ la with predictions from numerical simulations (Figure 
8l): iDuff v et al.1 120081: iGao et alll2012l: iBhattacharva e t al.l 

201a Diemer fc Kravtsov! 20141 : IMeneghetti et al.ll2014h . As 
noted above, the agreement between our results and the pre¬ 
dictions from 2013-2014 is excellent. More precise observa¬ 
tional constraints on the density profile shape of clusters, 
including mass dependence of the Einasto profile parame¬ 
ters await larger cluster samples, for example from the Dark 
Energy Survey (DES), Hyper Suprime-Cam survey (HSC) 
and the Large Synoptic Survey Telescope (LSST). 


120131 : iDiemer fc Kravtsov! 120141 ; IMeneghetti et ahl 120141 ') 
and individual clu ster mass measurements (Section 13.31) . 
lOkabe et al.l 1 20131 ') conducted a similar stacked lensing anal¬ 
ysis using a background galaxy catalogue ba s ed on a sin¬ 
gle colour cut, and based on the lllbert et al.l (120091 ') COS¬ 
MOS photometric redshift catalogue. The stacked shear sig - 
nal presented here is consistent with lOkabe et al.l (1201311 . 
The measurement uncertainties on the shear signal decrease 
as radius increases in this study due to radial dependence 
of our colour cut; this incre ases the weight of th e outer bins 
in our fit, relative to that of lOkabe et al.l (I2013T1 . Therefore, 
the mass and concentration from our new stacked analysis 
are marginally higher and lower than lOkabe et al.l d2013l l 
respectively. 

We compare the stacked result at A = 200, shown in 
Table [4] with the lognormal mean of the individual clus¬ 
ter mass and concentration measurements listed in Table [3j 
finding excellent agreement, with the latter being (M 200 ) = 
6.38 ± 0.24 x 10 14 /i _1 M 0 and (c 2 oo> = 3.73 ± 0.38. 

We also divide the clusters into two sub-samples of 25 
clusters based on the virial mass measured from the indi¬ 
vidual cluster shear profile models (Section 13.31) . adopting 
Mvi r = 8 x 10 14 as the dividing line between the 

sub-samples. We calculate the stacked shear profile for both 
sub-samples, and fit models, following the procedures ap¬ 
plied to the full sample. The results ar e in excell e nt agree¬ 
ment with both numerical simulations ( Bhattacharya et al.l 
120131 : IDiemer fc Kravtsovll2014l : IMeneghetti et al.ll2014l ~) and 

the best-fit mass-concentration relation for individual clus¬ 
ter mass measurements (Table [3] & |4] Figure [7]). 

_Some numerical simulation indicate that an lEinastol 

(119651 ') profile describes the spherically averaged mass den- 
sity profile for simulated halos better than the NFW profil e 
(I Navarro et al.l 120041 : IGao et all 120121 : iKlypin et al.l l2014h . 
The Einasto profile has the form: 


dlogp = _ 2 f r V 

dlogr \ r — 2 / 


(32) 


where r _2 is a scale radius at which the logarithmic slope 
is —2 and a is a shape parameter to describe the degree of 
curvature of the profile. The Einasto profile is specified by 
three parameters of M&, ca = r A /r _2 and a. We measure 
these three parameters by fitting the stacked lensing profile 
for all 50 clusters. As demonstrated by the NFW fitting, the 
contribution from the point mass is negligible compared to 
that of the main halo in the radial range 0.1 — 2.8 h -1 Mpc. 
We therefore just fit the Einasto profile and two-halo term. 
The best-fit parameters are a = 0.1611q q 4 i, M 200 = 
6.491^29 x 10 14 h~ 1 M@ and C 200 = 3.26lg;g4. These con- 


4 DISCUSSION 

In Section 14. II we quantify the remaining systematics in our 
analysis, in Section l4.2l we summarize our overall error bud¬ 
get, and in Section 14.31 we compare our mass measurements 
with results from the literature. 


4.1 Systematics 

In Section ld.l.ll we test the methods described in Section lXTI 
and that we use in Section 13. 2 1 to fit NFW models to the ob¬ 
served shear profiles. In Section 14.1.21 we correct the shear 
signal for the small colour selection and galaxy shape mea¬ 
surement biases calculated in Sections l2HlfelT7l and re-fit 
the NFW models to the corrected shear profiles. In Sec¬ 
tion 14.1.31 we calibrate the impact of using the full photo¬ 
metric redshift probability distribution of the background 
galaxies on our mass measurements. In Section 14.1.41 we 
consider the impact of forcing the number density profile 
of background galaxies to be flat. (We emphasize again that 
in our analysis and results we do not assume the number 
density profile to be flat.) 


4-1.1 Simulation Tests 


The radial range over which recent cluster weak- 


lensing studies (e.g. Ilsrael et al. 20121: 

Melchior et al.l 120141: 

lADDleeate et al. 

2014; Hoekstra et alJ 

2015f) have modeled 


the shear profile has been motivated in part by results from 
numer ical simulations. Here, we expand upon lOkabe et al.l 
J2013IL to test our individual mass measurements (Section 
El using synthetic weak shear catalogues based on simu¬ 
lated clusters and toy models. The former have the advan¬ 
tage of incorporating the full effects of the large-scale struc¬ 
ture that surrounds massive clusters, whilst the latter have 
the advantage of toy model clusters having perfectly known 
properties, and the properties of the background galaxy cat¬ 
alogues are matched to the observational data. Importantly, 
we calibrate the specific model fitting method that we ap¬ 
ply to our observational data directly on simulations whose 
properties match our own sample and data. 

We use mock observations of clusters from the “Cosmo- 
OWLS” cosmological hydrodynamical simulation that re¬ 
produces a large number of local gala xy cluster scaling re¬ 
lations, within a 400 /i -1 Mpc box d£e Brun et al.l 120141 : 
iMcCarthv et al.ll2014l ~). We use the simulations that include 
cooling, star formation, supernova feedback and AGN feed¬ 
back with a heating temperature ATheat = 10 8 K, known 
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M 200 [10 14 /zT 1 M 0 ] 


Figure 8. Mass-concentration relation (left) and ct-mass relation (right) for the Einasto profile obtained by the stacked lensing analysis. 
The cross denotes the best-fit parameters and the contours show the 68.3%, 95.4%, and 99. 7% confidence level s . Left : blue dashed, gree n 
dott ed and magenta d a shed- dotted lines denote the mass-concentration relation derived bv lDutton fe Maccidl j2014l l lDuffY _^t 2 jajj_l[200jl. 
and iMeneghetti et all J20l4), respectively. Right: blue dashed and green dotted lines are the a-mass relation from IPuttorTfc Macciol 
ll2014h and lGao^T*ahf LooT l, respectively. 


as the AGN 8.0 model. Weak-lensing catalogues compris¬ 
ing 100 galaxies per arcmin 2 were constructed following 
iBahe. McCarthy fc Kiri ul (|2012l ). Specifically, the mock ob¬ 
servations include the effect of shape noise, cluster substruc¬ 
ture and triaxiality, and correlated large-scale structure, 
and ignore uncorrelated large-scale structure and observa¬ 
tional effects such as uncertainties in galaxy shape mea¬ 
surements and redshifts. Note that ignoring uncorrelated 
large-scale structure is not expected to affect the measure¬ 
ment of possible biases in mass measurements v ia tests such 
as those described here JHoekstra et al.l 1201 ll ). The mass 
of the simulated clusters spans 5 x 10 14 /iy 3 1 Mg < M 200 < 
17 x 10 14 hjf M@ at a redshift of zi = 0.23 with the WMAP7 
cosmology of fl m ,o = 0.272 and 12 a = 0.728. We randomly 
extracted galaxies from the parent synthetic weak shear 
catalogues to match statistically the cluster-centric num¬ 
ber density profiles of colour-selected background galaxies 
in the observational analysis (Section l2.7ll . and fitted NFW 
models to the shear profiles following exactly the procedure 
laid down in Sections nnifcnoi This was repeated for 30 
realizations, each containing 23 simulated clusters. 

To quantify the mass measurement, we de¬ 
fine Ax in terms of the geom etric mean: Ax__= 
exp ((ln(A fit /Ainput))) — 1, foll owing lUmetsu et al, (2014), 
and similar to the methods of iBecker fc Kravtsov (2011). 
Here, Xnt is the best-fit mass or concentration. We recover 
the input M 200 and C200 from the numerical simulations 
with negligible bias. The mean bias on the mass and con¬ 
centration measurements across the full suite of realizations 
of the simulations is < 1 per cent. The scatter between mea¬ 
surements of the bias on mass using individual realizations 
is 4.8 per cent, which is comparable with the measurement 
uncertainty of 4.4 per cent on the bias from an individual 
realization. Likewise the realization-to-realization scatter 


in bias on concentration is 6.2 per cent, with a typical 
measurement uncertainty on an individual realization of 6.3 
per cent (upper three panels of Figure 0. 

We repeat this test using cluster density profile mod¬ 
els based on analytic NFW halos, and construct synthetic 
background galaxy catalogues that match the observed cat¬ 
alogues as closely as possible. For each of 50 analytic clus¬ 
ter profiles we adopt the observed positions of background 
galaxies and randomly draw a galaxy from the full back¬ 
ground galaxy sample across all 50 clusters, thus simulta¬ 
neously randomising the galaxy orientations, and matching 
statistically the source redshift distribution. The NFW pa¬ 
rameters are randomly chosen from the measured values for 
our cluster sample. We compute the synthetic shear pro¬ 
file for each of these 50 analytic clusters 10 times and fit 
an NFW model following the procedures laid down in Sec¬ 
tions nm & 13.21 Again, we recover the input masses and 
concentrations with < 1 per cent bias. The scatter between 
realizations is 4.3 and 8.9 per cent on mass and concentra¬ 
tion respectively, and the typical measurement uncertainty 
on individual realizations is 7.7 and 14.5 per cent on mass 
and concentration respectively (lowest panel of Figure [9j - 

In summary, we conclude that our shear profile fitting 
algorithm, a key feature of which is the adaptive choice of 
binning scheme, recovers the mean mass of our sample with 
negligible bias. 

4-1.2 Shear Calibration and Background Selection 

In Sections 12. 41 fc 17771 we developed methods to measure the 
shape of faint galaxies and select faint red galaxies as back¬ 
ground galaxies with small systematic biases of 3 and 1 per 
cent respectively; both acting in the sense that we slightly 
under-estimate cluster mass. Here we estimate how these 
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Figure 9. Calibration tests for weak-lensing mass measurements 
(M200;C200)- The top three panels and the bottom panel repre- 
sent the results based on the numerical simulations (diamonds 
iMcCarthv et all 1201 ll : iBahe. McCarthy fc Kin3l2012l . ; 23 clus¬ 
ters) and the toy models (circles; 50 clusters), respectively. The 
filled and open symbols denote the geometric means for the mass 
and the concentration, respectively. 


biases propagate through to the actually cluster mass mea¬ 
surements. 

The shape measurement bias is expressed as the multi¬ 
plicative shear calibration factor, m = —0.03, following the 
STEP programme. Given m, we therefore correct the mea¬ 
sured tangential shear signal by A£+ —> A£+(l + m)^ 1 and 
repeat the tangential shear fitting described in Section f3.II 
We express the comparison between the original masses (Ta¬ 
ble [2j) and the corrected masses we define Am in terms of 
the geometric mean: Am = exp ((hi(M° ris /MX orr )>) - 1). 
We find that the corrected masses are ~ 3 — 5% higher than 
the original masses, the range of values reflecting the non¬ 
linearity of tangential shear profiles (Tabic [5j). 

Turning to the colour selection of background galaxies, 
the right panel of Figure [3] shows that the contamination 
levels are 1 per cent at r 1.3 /i -1 Mpc and well below this 
level at r it 1.3 /i -1 Mpc. We therefore boost the shear sig¬ 
nals at r ^ 1.3 h -1 Mpc and re-derive the cluster masses, 
again following Section 13.11 Expressing the comparison in 
the same manner as above, we found that the masses cor¬ 
rected for contamination are within < 1 per cent of the 
original masses (Table 0. 

For completeness, we combine these two shear correc¬ 
tion terms in quadrature to give an effective multiplicative 
bias of m = —0.032, and evaluate the mass measurement 
bias of the combined shear calibration and contamination ef¬ 
fects. As expected, the bias is mainly attributed to the shear 
calibration, with the combined correction yielding results in¬ 
distinguishable from the pure shear calibration correction. 
Individual cluster masses based on the corrected tangential 
shear profile are given in Appendix [B] 


4-1.3 Mass Estimates with photometric redshift P(z) 

In Section f2.6l we adopted as the redshift of each faint galaxy 
in our sample, the median of the stacked posterior probabil¬ 
ity distribution of the nearest 100 neighbours in the (V—i)—i 
space of the COSMOS catalogue. We therefore essentially 
adopted a point estimate of the redshift of each of our galax¬ 
ies. However it is well known that the photometric redshift 
probability distribution of galaxies can be asymmetric, and 
present multiple peaks. The full photometric redshift prob¬ 
ability density function, P(z), fully describes such implicit 
systematic uncertainties. Indeed some recent studies have 
used the ful l P(z) for some clust ers in their weak-lensing 
sample (e.g. lApplegate et al.ll2014l ). 

Here, we test whether our method that ignores the full 
P(z) available from the COSMOS survey suffers any sig¬ 
nificant bias. In a similar vein to Section csa we estimate 
the full P(z) probability function of individual galaxies as 
an ensemble average of p COSMOS ( 2 ) for 100 neighbouring 
COSMOS galaxies in the colour-magnitude plane, 


p (*) = ]v-£^ 


COSMOS 


to- 


(33) 


Given the probability function, the tangential shear compo¬ 
nent can then be calculated as follows: 


(AE+XrO = 


^Ejn Sz[ ff+ I ™^cr,n(Zs)'lC n (z s )P(Zs)dz s 

En f z T w n(z s )P(z a ) dz s 


(34) 


The errors for the shape noise ( C g ) and the photometric 
redshift ( C s ) are estimated as: 


1 E„ f™ ^g,n^cT,n( z ^ w n( z s)P( z s )dz s 


C<> " 2 (E n f™Wn(Zs)P(z s )dz s ) 2 


and 


Cs,ij — 


E„ (g+,nZ C r,n(z s )w n (z s ) - (A T,+}) 2 P(z s )dz a 

E n !~P{z.)dz, 


respectively. 

We select the same background galaxies used in our 
main analysis and compute the tangential shear profiles us¬ 
ing Equation (1341) . The radial bins are chosen using the same 
method as in our main analysis. The difference in the tangen¬ 
tial shear components estimated by the single source redshift 
and P(z) is calculated using: 



((AS?) - (A < M )) 2 

(a?) a + K W ) 2 


(35) 


where <j\ = C g + C s . The deviation is, on average, A+ ~ 
0.05, with smaller deviations of A+ ~ 0.02 at small radii 
(r < 0.3 /toMpc) and larger deviations of A+ ~ 0.07 at 
large radii (r > 1 h -1 Mpc). This is because the errors at 
small radii are larger, due to a the relatively small number 
of background galaxies in bins of smaller solid angle. When 
we ignore the errors, the average deviations are still neg¬ 
ligible, A+ ~ 0.01 at r < 0.3 /i -1 Mpc and A+ ~ 0.002 
at r > 1 Ii _1 Mpc. We also compare the best-fit mass and 
concentration parameters (Figure 1101) . The two measure¬ 
ments are in excellent agreement, with geometric means of 
0.97 ± 0.05 and 0.96 ± 0.15 for M 200 and C 200 , respectively. 
We also made a background galaxy catalogue using the P(z) 
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Figure 10. A comparison of the mass (left) and the concentration (right) at A = 200, estimated by the single redshift, ( z s \ Section 13.211 
and the probability density function (P(z); Section |4.1.31 1. 


function and our radius-dependent colour-cut (Section 12.71) . 
and found that again the mass measurements do not change 
significantly (0.96±0.05 and 0.96±0.15 for M 200 and 0200). 
We conclude that with the current sample and data we are 
unable to detect any systematic difference between the mass 
and concentration measurements based on the mean of COS¬ 
MOS point estimates of the redshift of individual galaxies 
and the full COSMOS P(z) function. 


4-1-4 Boost factor 

Imperfect selection of background galaxies leads to under¬ 
estimated weak-shear signals because background galaxy 
samples are contaminated by faint cluster galaxies, the 
shapes of which present no lensing signal due to the clus¬ 
ter. The number density of cluster galaxies decreases with 
increasing projected cluster-centric radius. Thus at fixed se¬ 
lection method, the ratio of cluster galaxies to background 
galaxies, /mem, decreases with increasing projected cluster¬ 
centric radius, and thereby dilutes the shear signal more 
at smaller radii than at larger radii. Our approach is to 
vary the colour cut used to select background galaxies as 
a function of cluster-centric radius, simply exploiting the 
declining number density of cluster galaxies as a function 
of radius, witho ut invoking any physical assumptions. An 
alternative (e.g. iKneib et all 120031 : lApplegate et alJ 12014 
iHoekstra et alj 20151 ) is to correct the measured shear sig¬ 
nal by a factor 1 ± / m em, and thus to assume that the 
observed number density profile of background galaxies is 
flat. This correction is referred to as a boost factor. How¬ 
ever, the assumption of a flat observed number density pro- 
file of background galaxies ignores magnification bias (e .g. 
I Broadhurst. Taylor fe Peacockl 1 19951 : lUmetsu et al.l 12014 ) - 


i.e. the depletion or enhancement of the number density of 
background galaxies due to lensing magnification. 

We compare the boost-factor method with our methods 
that do not invoke a boost factor and instead rely on selec¬ 
tion of red galaxies to achieve ^ 1 per cent dilution across 
all radii. For this purpose, we construct background galaxy 
catalogues that suffer dilution by selecting red galaxies with 
a positive colour-offset from the red-sequence, AC > 0, 
and compute the stacked shear profile for this diluted sam¬ 
ple of galaxies. As expected, the amplitude of this shear 
profile is suppressed relative to the shear profile from our 
main analysis, with the suppression increasing to ~ 50% at 
~ 100 h _1 kpc from the cluster centre (left and central panels 
of Figure fill) . The suppression of the signal due to dilution 
appears to be negligible at ^ 1 /i -1 Mpc. The contaminat¬ 
ing population of faint cluster galaxies is seen clearly as an 
excess of galaxies at <j 1 h _1 Mpc in the stacked number 
density profile of galaxies selected as having AC > 0 (right 
panel of Figure fill) . In other words, the excess of number 
density profile is negligible beyond r2oo- The evidence indi¬ 
cates an internal consistency that the number density excess 
and the stacked-lensing mass estimate are consistent with 
each other. Note that the number density profile in Fig¬ 
ure im is calculated after masking the solid angle subtended 
by bright galaxies ( i' < 20) out to elliptical radii a factor of 
3 than the elliptical shape parameters of SExtractor - i.e. 
corresponding to the isophotal limit of detected objects. We 
also fully consider the finite field-of-view in the number den¬ 
sity calculation. To quantify the contaminating population 
we fit the function / mem = Aexp(— r/ro) to the measured 
number density profile (dashed blue curve in right panel of 
Figure [Ul) . We use this model to boost the measured shear 
signal by a factor 1 ± / mem (blue crosses in left and cen¬ 
tral panels of Figure fill) . It is clear that the boosted shear 
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Figure 11. Systematic mass measurement errors caused by the boost factor. Left: the stacked tangential shear profiles. Back diamonds 
and green squares denote the profile computed by our background catalogue and galaxies with AC > 0, respectively. Blue crosses and 
red circles are the profile corrected by the boost-factor and the boost-factor taking into account the magnification bias, respectively. 
Middle: the differential fraction of lensing signals from the signals for our background catalogue. The lensing signal corrected by the 
boost factor, with the assumption of the uniform background distribution, is significantly underestimated. Right:the normalized number 
density profile. Green squares shows a clear excess of the density profile of galaxies with AC > 0. The best-fit model is represented by 
the dashed blue lines. The dotted red line denotes the depleted number density profile due to the magnification bias. The solid black line 
is the di l ution model (Section IO- The dotted-dashed magenta line is the boost factor for MACS J0417.5—1154 from the WtG project 
llApplegatc et A~I|20 1 41 . showing a shallow distribution. 


signal underestimates the lensing signal that we detect by 
~ 40% on small scales, and > 10% at all radii interior to 
~ 1 fi _1 Mpc. Lens magnification is an obvious culprit for 
this apparent deficit of signal in the boost-factor-corrected 
shear profile. 

The galaxy-count is depleted by the magnification bias, 
as expressed by 

N(r;<m) = No(r;<m)ij,(r) 2 ' 53 ~ 1 , (36) 

where fjT 1 = (1 — k) 2 — \^\ 2 is the lensing magnification which 
expands the area of sky and enhances the flux of galax¬ 
ies, and s = dlog 10 No(m)/dm is a logarithmic count slope. 
Given the best-fit mass model derived by the stacked shear 
analysis (Table [4j), we calculate the number density profile 
expected from the magnification bias, assuming s = 0.15 
dUmetsu et al.l2014l l . This calculation shows clearly that the 
expected number density profile of background galaxies is 
not flat, showing a decline interior to ~ 1 /i -1 Mpc (dotted 
red curve in right panel of Figure HD- This indicates that 
the assumption of a flat background galaxy number density 
profile is incorrect, even on scales comparable with rsoo- 

Next, we boost the AC > 0 shear profile (blue squares 
in left and central panels of Figure HT1) by both the boost fac¬ 
tor, and the expected number density profile of background 
galaxies from the magnification bias calculation discussed 
in the proceeding two paragraphs. This boost factor and 
magnification bias corrected shear profile comes closer to 
recovering our measured shear profile, although it remains 
~ 10% lower interior to ~ 1 /i -1 Mpc (red points in left 
and central panels of Figure fill) . Clearly, the number den¬ 
sity profile based on an imperfect background selection is 
tightly coupled with the dilution effect and the magnifica¬ 
tion bias. It is therefore very difficult to break the degener¬ 
acy between the dilution effect and the magnification bias 
using the imperfect background catalogue. We also mention 


that the boost-factor gives rise to systematics in the source 
redshift because member galaxies in background catalogue 
have inadequate redshifts. Overall, our analysis in this Sec¬ 
tion indicates that application of a boost factor without con¬ 
sideration of lens magnification may cause systematic biases 
even on quite large scales up to ~ rsoo- 

4.2 Error budget 

We summarize the contributions to our overall error budget. 
The largest systematic bias in our measurements is caused 
by imperfections in our faint galaxy shape measurements, 
with a STEP-like multiplicative bias of m ~ —0.03 that is 
independent of galaxy size. We tuned our faint galaxy se¬ 
lection method to yield a contamination level of ^ 1 per 
cent, thus again giving a very small systematic bias. This 
low level of contamination renders the background galaxy 
samples essential pure, and thus not requiring correction. 
Nevertheless, for completeness, we correct the shear profiles 
for these two bias terms and obtain “corrected masses” (see 
Appendix) that are ~ 3 — 5 per cent higher than the “un¬ 
corrected masses” in Tabled In particular, after correction, 
our measurements of M500 increase by just 2.7 per cent on 
average. 

We have also investigated several other errors, none of 
which contribute to systematic biases. First, the typical un¬ 
certainty on the distance ratio Di a /D s for each individual 
background galaxy, computed from the COSMOS photo¬ 
metric redshift catalogue following the methods described 
in Section Ml is cr/3 ~ 13 — 28 per cent. These uncertainties 
are included in the covariance matrix when fitting the NFW 
model to the respective shear profiles fSection l3.ll) . Second, 
when exploring different shear profile binning schemes, and 
applying them to the observational shear measurements, we 
found that the M 500 values for individual clusters scatter by 
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Table 5. Systematic errors on the mean mass, (M&), for 50 clusters, in units of a percentage. 


Name 

expected 

M vir 

M 200 

M 500 

M 1000 

M 2500 

Shear calibration 

-3 

-4.1 

-4.2 

-2.7 

-4.8 

-5.3 

Colour selection 

-1 

0.2 

0.0 

-0.5 

-0.9 

-1.7 

ShearH-colour calibration 

-3.2 

-4.1 

-4.2 

-2.7 

-4.8 

-5.3 

Radial bins 

- 

±0.9 

±1.0 

±1.4 

±2.3 

±4.4 


±1.4 per cent. Third, when applying our shear profile bin¬ 
ning scheme to simulated data we recovered the true masses 
with 1 per cent bias and ~ 5 per cent scatter between dif¬ 
ferent realizations of the simulations (Section [TTTT[| . Fourth, 
we tested whether using a point estimator of the photomet¬ 
ric redshift of a galaxy in the COSMOS catalogue intro¬ 
duced any systematic bias with respect to using the full 
P(z) distribution. We found that the latter yields masses 
4 ± 6 per cent lower than the former. Given the poor preci¬ 
sion to which have been able to measure this possible bias, 
we ignore it here. However this deserves more detailed inves¬ 
tigation in future larger surveys. Whilst we have sought to 
minimize the assumptions in our analysis pipeline, we would 
like to draw attention to the dependence of our results on 
the COSMOS field. We have assumed that the galaxy popu¬ 
lation probed by observations of this field are representative 
of the universe as a whole. This assumption is relevant to 
the SHERA-based shape measurement tests described in Sec¬ 
tion 1231 and aspects of our analysis that rely on the COS¬ 
MOS photometric redshift distribution (Sections 12.61 fe PZTl) . 
Several many-filter and deep COSMOS-like calibration fields 
would be very beneficial for future deep and wide lensing 
surveys including LSST. 


4.3 Comparison of mass measurements with the 
literature 


Ongoing pro j ects, including LoCuSS, CCCP 


llHoekstra et ah _ 20151') . CLA SH llUmetsu et al.l 120141 ') 

and WtG ( Applegate et al] l2014h . conducted cluster 
weak-lensing analyses. Comparison of the cluster mass mea¬ 
surements between these different surveys is of paramount 
importance for cluster cosmology experiments. There are 
important differences between the approach taken by each 
survey to the cluster mass measurements. For example, 
CCCP and WtG implement the boost factor method 
(Section 14.1.41) to correct their shear signal. Another 
key difference is the prior adopted on the concentration 
parameter of the NFW halo model that is fitted to the 
data. For example, WtG fix the concentration parameter at 
C 200 = 4, whilst we allow it to be a free parameter in our fits. 
We summarize the key differences between the respective 
analysis methods in Table [6] and highlight key points in 
the following sections. Throughout these sections, when we 
compare cluster masses between two surveys we do so within 
the same radii so as to avoid errors caused by aperture 
mis-match. Also, all comparisons are done without applying 
the colour-selection and shear-calibration corrections to our 
shear measurements, discussed in Section [4.1.21 


4-3.1 LoCuSS - Okabe et al. (2010) 

lOkabe et al.l d2010al ') conducted weak-lensing analysis for 30 
clusters using Suprime-Cam data. As some of clusters lacked 
V — band data, we measured masses only for 22 clusters, 
defining as background galaxies those galaxies with colours 
are redder or bluer than those of the red-sequence. The red 
and blue colour cuts were chosen by eye based on the run of 
lensing signal with colour offset from the red sequence. We 
didn’t adopt any correction factor inherent in the shear cal¬ 
ibration and the profiles of the background number density. 
The error on the tangential shear measurements only took 
into account shape noise. We treated the mass and the con¬ 
centration as free parameters when fitting the NFW model 
to the tangential shear profiles. The masses derived from 
our new analysis presented in this article are 20% to 9% 
higher than obtained in 2010, with larger differences found 
for masses measured in larger radii (Table[6] Figure fl2l) . Our 
previous mass measurements were under-estimated due to 
issues in shear calibration and imperfect background galaxy 
selection. 


4-3.2 CCCP Hoekstra et al. (2015) 

The Canadian Clu ster Comparison Project (CCCP 
iHoekstra et al.l 120151 1 carried out weak-lensing mass mea¬ 
surement for 52 clusters using the Canada-France-Hawaii 
Telescope (CFHT). Their sample is defined by clusters with 
an ASCA temperature of ksTx > 5keV in the range of 
0.15 < 2 < 0.55. They used B— and R— bands for the 
first 20 clusters with the CFH12k camera and g' — and r' — 
bands for the other 32 clusters with the Megacam (here¬ 
after we refer to the redder band as r-band for both instru¬ 
ments). They used solely their r-band data to select galaxies 
at 22 < r < 25 as background galaxies. The number den¬ 
sity profile of these galaxies was found to increase towards 
the cluster centres. Hoekstra et al. modelled this excess of 
galaxies on a cluster-by-cluster basis, and boosted the lens¬ 
ing signals by a factor of l + / C ontam(r) assuming a flat back¬ 
ground number density profile of background galaxies, where 
/contam (r) is the fraction of contaminating galaxies obtained 
from their models. This boost factor increases their masses 
by 1 — 2 per cent. They adop t the mass-concentration rela¬ 
tion (iDutton fe Macciai2014l ') for their mass estimates. 

There are 21 clusters in common between the CCCP 
and LoCuSS samples. On average, the CCCP masses are 
~ 3% —6% lower than our masses (Table[6] Figure fl2l) . how¬ 
ever this difference is not statistically significant. We note 
that the radial range of their model fits is 0.5 — 2/iy 0 1 Mpc 
i.e. it includes scales on which we expect the magnifica¬ 
tion bias to affect the slope of the number density profile of 
background galaxies (Figure fill see also Ziparo et al. 2015). 
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Figure 12. Mass c omparisons of differ ent projects: the top-left (OklO; [OkabeetjrlJ 2010ajb top-right fCCCP iHoekstra et alJl2015h . 
bottom-left (CLASH I Umetsu et alj[2014h and bottom-right (WtG lApplegate et al.ll2Q14 b Masses are measured at the same radii in order 
to avoid the aperture-induced errors. The mass ratios for the overlap sample are computed by the geometric means. 


Nevertheless, we see good agreement between LoCuSS and 
CCCP mass measurements. 


4-3.3 CLASH - Umetsu et al. (2014) 

The Cluster Lensing And Supe rnova survey with Hub¬ 
ble (CLASH: lUmetsu et al.ll2014l ) conducted a joint shear- 
and-magnification weak-lensing analysis of a sample of 20 
galaxy clusters at 0.19 H z ^ 0.69, using imaging through 
multiple filters with Subaru/Suprime-Cam. They measure 
galaxy ellipticities using the KSB+ method, and calibrate 
the isotropic PSF correction for galaxies detected with high 
signal-to-noise ratio - i.e. using methods similar to our own. 
They also employ a correction factor 1/0.95 to account for 
residual shear calibration. Background galaxies are selected 


in a colour-colour pl ane, typically based on the BjRcz'- 
band filters, following iMedezinski et all d20ld h They do not 
employ a boost factor to compensate for contamination of 
their background galaxy catalogues. The halo concentration 
for the NFW model is treated as a free parameter. We com¬ 
pare the CLASH and LoCuSS masses for 5 clusters in com¬ 
mon, obtaining excellent agreement (Figure [l^] and TableJS]). 
Recently, Umetsu et al. (2015) have published joint strong 
plus weak-lensing mass measurements of the CLASH sam¬ 
ple. Their strong-|-weak-lensing masses are in a similar ex¬ 
cellent agreement with our measurements. 
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Table 6. Comparison with previous studies, M pre (< r A*" ) /M LoCuSS(< All masses are the NFW spherical masses measured within 

the overdensity radii defined by previous studies (OklO,CCCP, CLASH and WtG) or 1.5hff Mpc, in order to avoid the aperture- 
induced-errors in the mass estimates. We also summarize the mass-measurement method and approach to correction of biases for each 
project. “Method” denotes the tangential shear fitting ((}+) and a joint fitting using tangential shear profiles and the magnification bias 
( g.y & //). “Calibration factor” is the shear-calibration factor, with “Yes” indicating that such a factor was applied to the shear signal 
before fitting mass models, and “No” indicating otherwise. “Boost factor” is the correction factor by the number density caused by 
imperfect background selection - Yes/No indicates whether or not this factor was calculated and applied to the data, ca states whether 
the concentration parameter was a free parameter in the fit, or fixed, or scaling with the mass. “Radial bins” gives the choice of radial 
binning scheme for the fitting of the shear profile. “N” is the number of clusters in common between this study and each of the other 
studies. 


Name 

Method 

CalibrationBoost 
factor factor 

Radial 

bins 

CA 

N 

Virial 

200 

500 

2500 

l-.5h.7Q 1 Mpc 

This paper 

' 9+ 

No 

No 

Adaptive Free 

- 

— 

- 

- 

- 

— 

OklO 

9+ 

No 

No 

Fixed 

Free 

20 

0.81 ±0.07 

0.82 ±0.06 

0.85 ± 0.05 

0.91 ±0.08 

- 

CCCP 

9+ 

Yes 

Yes 

Fixed 

Scaling 

21 

0.94 ±0.06 

- 

0.95 ± 0.06 

0.97 ±0.07 

- 

CLASH 

g+ &£ g 

Yes 

No 

Fixed 

Free 

5 

0.99 ±0.14 

1.02 ±0.13 

1.07 ±0.12 

- 

1.03 ± 0.10 

WtG 

9+ 

Yes 

Yes 

Fixed 

Fixed 

17 

- 

1.14 ±0.09 

1.16 ± 0.12 

1.17 ±0.14 

1.15 ± 0.07 


4-3.4 WtG - Applegate et al. (2014) 

The Weighing th e Giants programme (WtG; 
lApplegate et al.l |201J) have conducted weak-lensing mass 
measurements for 51 X-ray luminous galaxy clusters in the 
redshift of 0.15 ;$ z kj 0.7, using data from Subaru/Suprime- 
Cam. They calibrated their shape measurement pipeline 
using the STEP simulations and used this calibration to 
correct their faint galaxy shape measurements down to 25th 
magnitude. Note that their shape measurements suffered 
significant noise bias at fainter magnitudes. We compare 
our masses with WtG masses based on colour-selected 
galaxies, and therefore concentrate further discussion of 
their methods on those measurements. 

WtG define two catalogues of background galaxies to 
compute tangential shear profiles. The first catalogue is 
used to calculate the contamination correction, namely the 
boost factor, wh i ch is shown by blue points in Figure 4 of 
lApplegate et al.l ( 2014 1. and the second catalogue is used 
to compute the tangential shear profile represented by or¬ 
ange points in the left panel of their Figure 4. The first 
catalogue is defined by employing magnitude and size cuts 
and exclusions of stars and galaxies that lie on the red- 
sequence of cluster galaxies. They found an apparent ex¬ 
cess in the number density profiles at small radii because of 
imperfect background selection. The excess is described by 

,/rriem (r) — Af me m/( Amerci “t“ (Vbkg) — /*500 CXp(l c/r-,00. X ) 
which is the ratio of the member galaxies divided by the to¬ 
tal number (member and background) of galaxies selected. 
Here, they assume a constant number density profile of 
background galaxies, and thus ignore the potential effects 
of magnification bias. All clusters are fitted simultaneously 
to estimate /spp and rsoo.x- The normaliz ation /500 for the 
first c atalogue is (8.6 ± 0.9)%. Figure 4 in lApplegate et ahl 
(1201411 shows that rspp.x ~ 1.4 — 2/t 7fl 1 Mpc which is very 
large. For example. iMartino et al.l (12014 ~) obtain a median of 
r 5 oo,x — I hff Mpc and values in the range 0.8 — 1 ,6hff Mpc 
for the same sample that we study here, that has consider¬ 
able overlap with the WtG sample. Whilst cluster-by-cluster 
comparison of r$oo,x measurements is required to be cer¬ 
tain, the large values of this parameter used by WtG sug¬ 


gests that they assume a rather shallow radial density distri¬ 
bution when they model contamination of the background 
galaxy samples (right panel of Figure If 11) . Their shallow 
number density profile indicates that cluster members are 
distributed beyond the cluster virial radii, which conflicts 
with our results that the number density of cluster members 
is negligible in these regions (Section 14. 1.41) . Their number 
density result therefore appears to be inconsistent with their 
weak-lensing mass measurements. The shallow number den¬ 
sity profile found by the WtG may therefore be a source of 
systematic bias in their mass measurements. 

Next, they draw the second catalogue from the first cat¬ 
alogue with additional conditions of robust shape measure¬ 
ments and lensing cut ( S/N > 3 and > 1.15rh*)- The tan¬ 
gential shear profile is computed from the second catalogue. 
The excess of the number density profile, ft ,pp = 4.8 ± 1.6%, 
is less than that obtained by the first catalogue. The tangen¬ 
tial shear profile using the second catalogue is corrected by 
the boost factor that was calculated using the first catalogue. 
We here explicitly describe their definition of the corrected 
lensing signal, as follows, g+ nd —> ff+ nd /(l - /mem(r-)). Here, 
1st and 2nd denotes the quantities computed by the first 
and second catalogues, respectively. This inconsistency may 
affect the precision of the WtG mass measurements. 


We emphasize that the boost factor cannot recover the 
correct lensing signal because of the magnification bias, 
as demonstrated in Sections vm and 14.1.41 notwithstand¬ 
ing the fitting range, 0.75 — -ihff Mpc, that WtG adopt. 
Finally, they assumed C2pp = 4 for the NFW model for 
all clusters at all redshifts. Numerical simulations (e.g. 
Buflocke^jd. 2001 j_ DuffvjTaL [200^_ Bhattachaj^et_al] 


20131 : iDiemer fe Kravtsov! 120141 ; iMeneghetti et al.l 120141 : 
Ludlow c t al.l 20141 1 show that the concentration depends 
on both the halo mass and redshift. 


WtG masses are ~ 15% higher than our masses, in¬ 
dependent of overdensity, albeit at ~ 1 — 2cr significance 
(Tabic [6] Figure fl2ll . Note that we exclude A1758N from 
this comparison because the WtG adopt a radial fit range 
that extends into the companion cluster A1758S. When we 
follow the WtG method - i.e. restrict the radial range of 
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Figure 13. Mass comparisons for 12 cluster over lapped 
among the three proje cts of CCCP llHoekstra et al .1 \‘20 1 5li . WtG 
IIApplcgatc et ajj I2014T1 and LoCuSS. Masses are measured at 
1 ,5hy 0 l Mpc, derived by the WtG fitting method; the radial range 
is 0.75 — 3 hyy Mpc and C 200 = 4. The bottom panel shows a 
comparison of CCCP masses (red diamonds) and WtG masses 
(blue circles). The top panel shows individual masses of 12 clus¬ 
ters. Red diamonds, green squares and blue circles denote CCCP, 
LoCuSS and WtG masses, respectively. 


the fit to 0.75 — Sh^Mpc and adopt C200 = 4, the statisti¬ 
cal significance of the disagreement increases slightly, with 
geometric means at A = 200,500,2500 and 1.5 hf,/Mpc of 
1.18 ± 0.09, 1.19 ± 0.12,1.19 ± 0.14, and 1.18 ± 0.07, respec¬ 
tively. We note that the difference between our mass mea¬ 
surements and those of WtG is significant at > 2a on the 
scale preferred by WtG, 1.5 Mpc. 

In summary, we have identified several strong assump¬ 
tions and inconsistencies in the WtG analysis: the incon¬ 
sistent calculation for the boost factor, the shallow num¬ 
ber density profile leading to high boost factor, ignoring the 
magnification bias, and the fixed concentration C200 = 4. We 
expect that the tension between our respective mass mea¬ 
surements will be caused by one or more of these issues. 


4.3.5 CCCP, WtG and LoCuSS 

The cleanest comparison between surveys is between the 
respective mass measuremen ts for 12 clus t ers i n com¬ 
mon between all of CCCP dHoekstra et al.l 120151) . WtG 
(I Applegate et al.1l2014l ) and LoCuSS (this article). We fol¬ 
low the WtG fitting method, that is, we fix the NFW model 


concentration parameter at C200 = 4, and fit the NFW 
model to the observed shear profile in the radial range is 
0.75 — S/i^Mpc. We compare three spherical NFW masses 
within l.S/i^Mpc. Based on these conditions, the mass com¬ 
parison approximately corresponds to a comparison of lens- 
ing signal at the fixed radial range. This like-for-like com¬ 
parison confirms the results discussed in the preceding sec¬ 
tions: for the subsamples of clusters in the LoCuSS sample 
that have also been studied by CCCP and WtG, the CCCP 
and LoCuSS masses are in good agreements with each other, 
while the WtG masses are ~ 2a higher than both CCCP and 
LoCuSS masses (Figure fl3l) . 


4-3.6 Sensitivity of conclusions to “sample selection” 


An important caveat on the conclusions discussed in Sec¬ 
tions [43L2JJ43L5] is that they are strictly only applicable to 
the specific sub-samples of clusters that have been observed 
by other surveys that are in common with our sample. We 
therefore consider whether our conclusions are supported by 
inter-survey comparisons in the literature - i.e. by compari¬ 
son of masses of clusters in the overlap between other pairs 
of surveys, thus excluding our LoCuSS sample selection and 
mas s measurements f rom t he discussion. 

_ Hoekstra et al. d2015l ) find that WtG masses from 

lApplegate et al.l ( 20141) are ~ 6 — 8 per cent higher than 
their CCCP measurements at ~ 2a significance, depend¬ 
ing on the redshift distribution adopted by the two surveys. 
This is consistent with our conclusions, however it is intrigu¬ 
ing to note that WtG masses agree very well with CCCP 
masses when the lat ter are compute d by deprojecting aper¬ 
ture mass estimates. iHoekstra et ahl also compare their mass 
measurements with CLASH measurements of 6 clusters in 
com mon between these su rveys, finding that CLASH masses 
from lUmetsu et alj (120141) exceed CCCP masses by 12 ± 5 
per cent. Umetsu et al. (2015) also find similar results (a 
16 ± 10 per cent excess) for five clusters in common between 
their joint stron g plus weak-lens ing analysis of the CLASH 
sample and the Hoekstra et al.l weak-lensing masses. The 
CCCP/CLASH comparison appears to be at odds with our 
finding that CLASH, CCCP, and LoCuSS are all consis¬ 
tent, and underlines our caveat that conclusions based on 
mass comparison between surveys may be sensitive to the 
specific overlapping sub-samples considered. However, the 
overlap between CLASH and WtG numbers 17 clusters, i.e. a 
number of clusters comparable with the LoCuSS/CCCP and 
LoCu SS/WtG samp les discussed in Sections 14.3.21 & 14.3.41 
Both lUmetsu et alj J2014I ) and Umetsu et al. (2015) find 
that WtG masses exceed the CLASH masses by ~ 7 — 10 
per cent, albeit at ~ la significance. 

In summary, a consistent picture emerges from compar¬ 
isons between LoCuSS, CLASH, CCCP, and WtG surveys 
when the samples in common between the surveys number 
> 10 objects: the WtG cluster mass calibration exceeds that 
of other surveys by ~ 6 — 15 per cent at ~ 1 — 2a significance 
per pair of surveys considered. Further inter-comparison of 
mass measurements based on all of the methods used by the 
respective surveys, preferably using an enlarged and com¬ 
mon sample of clusters, should help to clarify the remaining 
differences between the surveys. 
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5 SUMMARY 

We observed an almost complete Lx-selected sample of 50 
galaxy clusters at 0.15 < 0 < 0.3 through V- and i'-band fil¬ 
ters with the Suprime-CAM instrument on the Subaru 8.2-m 
telescope. We used these data to measure the weak gravi¬ 
tational shear signal, and thus to infer the total mass and 
concentration of each cluster, and the sample as a whole via 
a stacking analysis. The size of our sample and typical sta¬ 
tistical precision on weak-lensing mass measurements of 30 
per cent motivates our goal of controlling systematic biases 
in our analysis at the ~ 30/V50 ~ 4 per cent level. 

The recent literature identifies the dominant systematic 
uncertainty in weak-lensing analysis of cluster samples as 
contamination of background galaxy samples by faint cluster 
galaxies. We extend our background galaxy selection method 
based on the dependen ce of gravitationa l shear signal on 
the colours of galaxies (lOkabe et alj 12013H to incorporate 
the empirical fact that, at fixed colour cut, contamination 
by faint cluster members is a declining function of cluster¬ 
centric radius. This allows us to define colour cuts that are 
a function of clustercentric radius such that we achieve neg¬ 
ligible contamination (we required ^ 1 per cent contamina¬ 
tion) whilst achieving a number density of 13 galaxies per 
square arcminute. We stress that our approach to selecting 
background galaxies neither assumes a form for the clus¬ 
ter mass distribution nor assumes a radial distribution for 
the background galaxies. We show that this latter point is 
important because gravitational magnification modifies the 
observed radial distribution of background galaxies on scales 
as large as > 1 /i _1 Mpc from the cluster centres and thus 
complicates methods that assume the radial distribution of 
background galaxies is flat (see also Ziparo et al. 2015 for a 
detailed discussion). 

The dominant systematic bias in our analysis is the ac¬ 
curacy of our faint galaxy shape measurements. We test our 
modified KSB+ method on large simulated datasets that 
incorporate smooth analytic galaxy templates based on Ser- 
sic profiles, and realistic galaxies based on high resolution 
imaging with the Hubble Space Telescope. The key feature 
of our shape measurement code is that we use a very high 
signal-to-noise ratio cut when selecting galaxies to model 
the isotropic PSF correction. The precision that we achieve 
in this model as a result of this cut is important because it 
significantly reduces the noise bias t hat is typically seen at 
faint flux levels in ot her studies (e.g. lApplegate et al.ll20l4 
iHoekstra et al .1 [2015)1 ). We show that our STEP-like multi¬ 
plicative shape measurement bias is of order 3 per cent, and 
is not a strong function of apparent magnitude and galaxy 
size. This is the largest systematic bias that we have identi¬ 
fied in our analysis. 

Controlling contamination of our background galaxy 
sample at the per cent level across a wide range of cluster¬ 
centric radii (down to ~ 200 /i _1 kpc) affords us the opportu¬ 
nity to fit for both mass and concentration parameter when 
modeling the weak shear signal. This is important because 
it eliminates the possibility of systematic biases caused by 
fixing the shape of the density profile or adopting a mass con¬ 
centration relation from simulations. We are also interested 
in measuring the mass-concentration relation of cluster-scale 
dark matter halos. To further mitigate the possibility of bi¬ 
ases when modeling the shear signal, we apply a range of 


binning schemes (number of bins, inner and outer fit radii) 
to the data that are well motivated by the physical prop¬ 
erties of clusters and previous tests of weak-lensing meth¬ 
ods (e.g . IMeneghetti et akll20ld : iBecker fe Kravtsovl l201ll : 
iBahe. McCarthy fc KiruJ 2012h . For each cluster we adopt 
the mass and concentration measurement of the binning 
scheme that yields the measurements that are the closest 
to the mean of the measurements from the range of binning 
schemes that we explore. This is a new method, and thus we 
test it on simulations including full hy drodynamical numer¬ 
ical simulations from cosmo-OWLS jLe Brun et alj 120141 ; 
iMcCarthv et al.l 120141 4 , based on an AGN model that re¬ 
produces a large number of local X-ray-SZ-optical scaling 
relations. Our tests show that we recover the ensemble mass 
and concentration of the simulated systems with sub per 
cent accuracy and of order 5 per cent scatter between differ¬ 
ent realisations of the simulated lensing observations. Issues 
relating to the binning and modeling of the shear signal are 
therefore sub-dominant in our analysis. 

We now summarise our main science results. We mea¬ 
sured weak-lensing mass for individual clusters by fitting 
tangential shear profile with the spherical NFW model and 
investigated the mass-concentration relation. The best-fit 
mass-concentration relation is i n excellent agreemen t with 
recent numerical simula tions (|Bhat4jMiaj^i^L_alJ 120131 : 
iDiemer fe Kravtsovll20l3 : IMeneghetti et alT 20l4) . We also 
measured the average mass density profile for the NFW and 
Einasto profiles by stacked lensing analysis, considering the 
point mass and the two-halo term. The best-fit NFW model 
of the stacked signal agrees well with the mean of the in¬ 
dividual mass measurements, supporting that the stacked 
lensing analysis recover the average of individual mass mea¬ 
surements. The Einasto profile also agrees with numerical 
simulations, albeit large scatter between numerical simula¬ 
tions. 

We compared our lensing masses with masses from 
other projects (CCCP, CLASH and WtG). The philos¬ 
ophy of weak-lensing analyses for the four projects in¬ 
clude some strong difference s. Our mass measure ments agree 
within lcr with the CLASH dUmetsu et al.ll2014T ) and CCCP 
dHoekstra et al.|[2015l f surveys, with whom we have 5 and 21 
clusters in common respectively. Our mass measurements 
are ~ 15 per cent lower than WtG measurements for the 
17 clusters in common, at ~ 1 — 2 a significance, depending 
on the mass measurement aperture. The fairest comparison 
between CCCP, WtG and LoCuSS is for the 12 clusters in 
common between all three surveys and based on a common 
mass measurement aperture and assumption on density pro¬ 
file shape. Adopting C200 = 4 and measuring mass within 
1.5hjf Mpc (i.e. matching to WtG), confirms that LoCuSS 
and CCCP mass calibrations are consistent with each other 
and WtG is in tension at 2 a higher mass than the other two 
surveys. 


To guard against possible sensitivity of this conclusion 
to the specific samples under investigation here, we also re¬ 
viewed recent studies in the literature that compare other 
pairs of surveys (and thus exclude LoCuSS). We find that 
WtG are in tension with masses ~ 6 — 15 per cent higher 
than other surveys at ~ 1 — 2 a significance in all pairwise 
comparisons that include )> 10 objects. We discuss a range 
of possible causes of this tension, including strong assump¬ 
tions on the radial distribution of background galaxies and 
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the halo concentration parameter, and several inconsisten¬ 
cies in the construction of the WtG background galaxy cat¬ 
alogues. We expect that the tension is likely caused by one 
or more of these factors. 

In summary, we have controlled systematic biases in 
our weak-lensing analysis of a large sample of clusters at 
0.15 ^ a < 0.3 at the 4 per cent level. Therefore, as far as we 
can tell, systematics do not dominate our results. Our meth¬ 
ods have numerous innovative features that set us apart from 
contemporary surveys. In that context it is encouraging that 
out of four surveys including our own, three agree within la 
and one is in tension with the other three at just ~ 1 — 2cr. 
This represents important progress towards convergence on 
the mass calibration of galax y clusters for cosm ological sur¬ 
veys. In a companion article JSmith et al IB) we compare 
our weak-lensing mass measurements with estimates based 
on the assumption that the intracluster medium is in hy¬ 
drostatic equilibrium, and discuss the recently reported ten¬ 
sion between the Planck results from the primary CMB and 
galaxy cluster counts. 

Looking to the future, we consider the development of 
additional deep, i.e. > 26th magnitude, COSMOS-like pho¬ 
tometric calibration fields to be vital for further progress 
on cluster mass calibration. This will be particularly impor¬ 
tant to make secure progress on the calibration of clusters at 
higher redshifts, especially given that upcoming surveys con¬ 
centrate on a relatively small number (typically 4-6) of pho¬ 
tometric filters. Further testing of faint galaxy shape mea¬ 
surement techniques in the high shear regime and down to 
faint photometric limits, preferrably in collaboration with 
experts in the cosmic shear community, will also be very 
helpful. Overall, the emerging consensus between surveys at 
2 < 0.3 that use very different methods, encourages us that 
the future for cluster cosmology is bright, and we look for¬ 
ward to further progress from ongoing/future optical/near- 
infrared surveys including KIDS, DES, HSC, Euclid, and 
LSST. 
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APPENDIX A: RADIAL POSITION FOR 
TANGENTIAL SHEAR PROFILE 


We consider how to compute the mean tangential shear av¬ 
eraged over background galaxies located in the i-th bin of 
a shear profile, spanning the radial range n < n < r^- It 
is important to compute these bin radii accurately so at to 
minimize systematic biases in cluster mass measurements. In 
the limit of a uniform sheet of background galaxies with no 
intrinsic ellipticity and an infinite number density of back¬ 
ground galaxies, the average tangential shear in the i-th bin 
can be calculated analytically in the continuum limit: 


. Zl g+(r)rdr 

{9+)i ~ prdr 

J 7*1 


(Al) 


If one assumes a singular isothermal sphere (SIS) model, 
g+ ~ 7 + = Jr -1 , then in the weak limit, the average tan¬ 
gential shear in i-th bin becomes: 

(g+h = ( A2 ) 


where r m jd = (ri + r 2)/2 - i.e. the average tangential shear 
corresponds exactly with that of a SIS at a radius that is 
the mean of the inner and outer radius of the bin. For com¬ 
pleteness, we also use the SIS model to demonstrate that 
the correct way to calculate the radius at which to place a 
binned shear measurement is to weight by the shear signal, 
thus: 


Ti = 


fri 9+ r dr 


(A3) 


which yields rt = r m id for the SIS model, as above. So far 
we have used the SIS model as a well-motivated illustra¬ 
tion in the sense that tangential shear profiles of clusters 
do not deviate significantly from a power law slope of r . 
It is therefore interesting to note the convenience of using 
T'mid because it is trivial to compute and is independent of 
the details of the trial mass distribution when computing 
the likelihood of a model, under the assumptions described 
above. 

However the observed density profiles of galaxy clusters 
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Figure Al. The deviation between the averaged lensing strength 
and the input tangential shear measured at three different radial 
definitions. The black diamonds and while circles are the average 
over all radial bins and > IMpc, respectively. The top and lower 
middle panels show large and small numbers of background galax¬ 
ies in the large FoV. The upper middle panel is based on a large 
number of background galaxies with the Suprime-Cam’s FoV, cor¬ 
responding to stacked lensing analysis. The bottom panel is based 
on a small number of background galaxies with the Suprime- 
Cam’s FoV, corresponding to individual cluster lensing analysis. 


do deviate from the SIS model (e.g. lOkabe et al.l 120131 ). It 
is therefore important to check the expected amplitude of 
systematic biases incurred if r m id is adopted for clusters that 
do not have a SIS density profile. We therefore consider a 
power-law profile, p+ = Ar~ a , in which case integrating 
Equation IaTI yields an average tangential shear in a radial 
bin of: 


(g+)i 


(2 - a) (rl - r?) ' 


(A4) 


The slope, a, for the NFW model is equivalent to ~ 0.5 on 
small scales and 1.5 on large scales. To estimate the possible 
systematic bias caused by adopting r m id as the radius of a 
binned shear measurement, we define n as: 


A = 


( 2 -a)(r t-rl) 
2(r 2 2 -“ - r?-“) ’ 


(A5) 


and calculate the deviation of Vi from r m id for the cases 
a = 0.5 and a = 1.5. We find that the deviation of ry from 
r m id is less than 1 per cent in these cases. We therefore con¬ 
clude that, in the idealised case of uniform distribution of 
background galaxies, infinite number density of background 
galaxies, and no intrinsic ellipticity, r m jd is a sufficiently ac¬ 
curate radius for binned shear measurements, if the goal is 
to control systematic biases at the per cent level. 

In reality, background galaxies are intrinsically ellipti¬ 
cal, have a finite number density, and non-uniform distribu¬ 
tion on the sky. We therefore now consider how to relate the 
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idealized case to a real observation. We begin by discretizing 
Equations 1 All & IA3I 


(g+)i 


En g+,n w n 

E„«± 


(A 6 ) 


En 9+,n 

E ra 9+,n W n 


(A7) 


Note, when comparing Equations IA6I fc 171 that here we as¬ 
sume for simplicity that all background galaxies are at the 
same redshift. The previous discussion of the idealized case 
encourages the adoption of g+ oc A 1 in the realistic case, 
which yields a definition of r, as the weighted harmonic 
mean radius of the background galaxies in a given bin: 


f whm 


XX W n 


(A 8 ) 


This definition has the important advantage that it does not 
rely on knowledge of the cluster density profile in order to 
compute the shear profile, and is therefore simple and sta¬ 
ble to implement. For example, one alternative would be to 
use the measured g+ for each galaxy in Equations I A 6 I fc I A7I 
However the intrinsic ellipticity of galaxies results in this ap¬ 
proach producing very noisy results, in some cases returning 
r; < 0. A second alternative is to recompute the shear profile 
at every step when fitting a model to the shear profile, with 
the radii of the bins adjusted to reflect the shape and am¬ 
plitude of the density profile at each step. This approach is 
computationally expensive, and given the level of accuracy 
required, it is not currently justified. 

We test three alternative definitions of rt for tangen¬ 
tial shear profiles in order to validate the arguments set out 
above for a realistic case. We conduct a Monte-Carlo simu¬ 
lation using a mock catalogue of 1000 clusters. We assume 
a spherical NFW mass model. The cluster masses are uni¬ 
formly and randomly distributed in the log in the range 
2 x 10 14 M 200 < 2 x 10 15 h~ 1 MQ. The concentra¬ 

tion of each cluster is calculated fr om the respective cluste r 
mass via the mass-concentration (IMeneghetti et aDl20l4l . 
taking into account intrinsic scatter. The cluster redshifts 
are uniformly and randomly distributed in 0.15 < z < 0.3. 
The background galaxy redshifts are fixed z 3 = 1. We fixed 
the number of radial bins IVbin = 6 for the tangential shear 
profile in the range 100 — 2800 /i^q 1 kpc. We calculate the de¬ 
viation between the average measured tangential shear and 
the input tangential shear, measured at r m id and r w h m , as 
defined above. We also employ a third definition of radius 
for comparison, the area-weighted radius, defined as follows: 


fwgt 


En W ^ r n 

T, n Wn 


(A9) 


that has been used by several groups (e.g. iNiikura et al.l 
l2015l l. The deviation is averaged over all radial bins and 
over those bins at r > 1 hjf Mpc. 

The Monte Carlo simulations consider four observing 
scenarios. First, a large field-of-view (FoV; 60' x 60') and 
large number of background galaxies (nbkg = 300 arcmin -2 ), 
which is approximately equivalent to the idealized case dis¬ 
cussed above (top panel of Figure lATl) . Second is the FoV 
of the Suprime-Cam (24' x 37 / ) and nbkg = 300 arcmin -2 , 
corresponding to the stacked lensing analysis (Section 13.41 
upper middle panel of Figure 1X1)1 . Third is the large FoV 


and the small number of background galaxies (nbkg = 
10 arcmin -2 ); lower middle panel of Figure lATl Fourth is the 
FoV of the Suprime-cam and nbkg = 10 arcmin -2 , represent¬ 
ing our analysis of individual clusters (Section 13.21 bottom 
panel of Figure [ATI 

As expected, based on the discussion above, in the ideal 
case (top panel) case the lensing signal measured at both 
r m id and r w h m recover the input values accurately, while the 
signal measured at r W gt is biased high by ~ 1 per cent. In 
the case that represents our stacked weak-lensing analysis 
(upper middle panel) the shear measured at r m id suffers ~ 9 
per cent bias on large scales due to the finite FoV impinging 
on the outer annuli. On the other hand, the performance of 
v whm and r wg t are invariant to the size of the FoV, draw¬ 
ing attention to how these radii automatically take account 
of departures from azimuthal symmetry in the bin defini¬ 
tions. Similar results are seen in the bottom panel. Overall, 
the weighted harmonic mean radius (r w h m ) well describes 
the average tangential shear signal in all setup conditions, 
while the area-weighted radius r wgt gives biases of a few per 
cent in all cases, and r m id is vulnerable to the accuracy of 
corrections for annuli that are not fully covered by the ob¬ 
servational data. We therefore conclude that the weighted 
harmonic mean radius is the lowest bias radius at which to 
place the binned tangential shear measurements for the tan¬ 
gential shear profiles in individual and stacked weak-lensing 
analyses. 


APPENDIX B: MASS ESTIMATES WITH 
CORRECTION FACTORS 

Cluster mass calibration is of primary importance for 
cluster-based cosmology. Although we do not apply any cor¬ 
rection factor in our analysis, the systematic uncertainty 
(m = —0.03; Section [TTT 2 ]) inherent in the shear calibration 
and the 1 per cent contamination in our background galaxy 
selection might be not negligible for some scientific motiva¬ 
tions like precision cosmology. We therefore tabulate clus¬ 
ter masses determined by fitting the tangential shear profile 
corr ected with th e syste matic uncertainty in Table iBlI Note 
that ISmith et all d 201 6ll based their analysis on the masses 
listed in the submitted version of this p aper that are slightly 
different from those listed in Table [BT1 Smit h et aJJs result 
of /3p = 0.95 ±0.04 and /3x = 0.95 ±0.05 is unchanged when 
the updated masses in Table ED are used. 

This paper has been typeset from a TjgjX/ DT)bX file prepared 
by the author. 
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Table Bl. Mass measures for individual clusters on (M^) CO g determined by lensing signals with corrections of the shear calibration and 
the contamination of member galaxies; g+/( 1 +- m) where m = —3.2%. See the details in Section 14.1.21 
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